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Abstract 



Although in the case of polymer solutions the existence of hydrodynamic screening has been theoret- 
ically established some time ago, use of the same methods for suspensions of hard spheres thus far have 
failed to produce similar results. In this work we reconsider this problem. Using superposition of topologi- 
cal and London-style qualitative arguments we prove the existence of screening in hard sphere suspensions. 
Even though some of these arguments were employed initially for treatments of superconductivity and 
superfluidity, we find analogs of these phenomena in nontraditional settings such as in colloidal suspen- 
sions, turbulence, magnetohydrodynamics, etc. In particular, in suspensions we demonstrate that the 
hydrodynamic screening is an exact analog of Meissner effect in superconductors. The extent of screening 
depends on the volume fraction of hard spheres. The zero volume fraction limit corresponds to the normal 
state. The case of finite volume fractions-to the mixed state typical for superconductors of the second 
kind with such a state becoming fully "superconducting" at the critical volume fraction (p* for which the 
(zero frequency) relative viscosity '^{relative) diverges. Brady and, independently, Bicerano et al using 
scaling-type arguments predicted that for tp close to tp* the viscosity r]{relative) behaves as C{l — ip/(p*)~'^ 
with C being some constant. Their prediction is well supported by experimental data. In this work we 
explain such a behavior of viscosity in terms of a topological-type transition which mathematically can be 
made isomorphic to the more familiar Bose-Einstein condensation transition. Because of this, the results 
and methods of this work are not limited to suspensions. In the concluding section we describe other 
applications ranging from turbulence and magnetohydrodynamics to high temperature superconductors 
and QCD, etc. 
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1 Introduction 



In his 1905-1906 papers on Brownian motion for suspensions of hard spheres, Einstein obtained now fa- 
mous relation for the self-diffusion coefficient Dq for the noninteracting hard spheres of radius R immersed 
in a solvent at temperature T [1]: 

Do='-^ = ^. (1.1) 
7 6wRriQ 

In this formula Tyg is the viscosity of a pure solvent and fc^ is the Boltzmann's constant. This result is 
valid only in the infinite dilution limit. In another paper [2] Einstein took into account the effects of finite 
concentrations and obtained the first nonvanishing correction to tiq for small but finite concentrations. It 
is given by 

r?/r?o = 1 + 2.5 + O (^2) (1.2) 

with if being the volume fraction := y^nR^. In this formula n is the number of monodisperse hard 
spheres in the volume V. If wc formally replace rjQ by 77 in Eq. (1.1), the obtained result can be cautiously 
used as a definition for the cooperative diffusion coefficient D, i.e. 

D = ^. (1.3) 

Below, we use symbols Dq for the self- diffusion coefficient and D for the cooperative diffusion coefficient. 
By combining Eq.s (1.1) - (1.3), we also obtain: 

D/Do = l-2.5ip + 0{ip'^). (1.4) 

Eq.(1.4) compares well with experimental results, e.g. those discussed in Ref.[3]^. Numerous attempts 

have been made to obtain results like Eq.(1.4) systematically. The above results arc restricted by the 
observation that Stoke's formula for friction 7 is applicable only for time scales longer than the charac- 
teristic relaxation time r^. of the solvent, r^. := pR^/nrjQ), e.g. see [4]. In this formula p is the density of 
pure solvent. This requirement provides the typical cut-off time scale, while the parameter R serves as a 
typical space cut-off for the problems we are going to study in this work. 

By analogy with the theory of nonideal gases, expansion Eq. (1.2) is referred to as a "virial". Unlike 
the theory of nonideal gases, where the virial coefficients are known exactly to a very high order [5], 
values for coefficients in the virial expansion for rj have been an active area of research to date even in the 
low concentration regime. A considerable progress was made in obtaining closed form approximations 
describing the rheological properties of suspensions of hard spheres in a broad range of concentrations 
[6-8]. Similar results for particles of other geometries are much less complete [7,9]. An extension of these 
results to solutions of polymers has taken place in parallel with these developments [10]. A noticeable ad- 
vancements have been made in our understanding of rheology of dilute and semidilute polymer solutions 
for fully flexible polymers and rigid rods. It should be noted, though, that polymers add further com- 
plexities because the connectivity of the polymer chain backbone plays an essential role in calculations of 
rheological properties of polymer solutions. The effect of chain connectivity on viscoelastic properties of 
polymer solutions has been an object of extensive discussion, and many of theoretical difficulties encoun- 
tered in describing these solutions are shared by suspensions of hard spheres. In particular, it is known 
[11], that particles immersed in a viscous fluid affect the motion of each other both hydrodynamically and 
by direct interaction (hard core, etc). Since the motion of particles in a fluid is correlated, it contributes 
to the distribution of local velocities within the fluid. Behavior of many systems (e.g. those listed in 
Section 6) other than the hard sphere suspensions happens to be closely related or even isomorphic to 
that noticed in suspensions. This observation makes study of suspensions important in many areas of 
physics, chemistry and biology. For reasons which will become apparent upon reading, in this work we 
shall mention only physical applications. 

In a polymer solutions when the polymer concentration (p increases, it is believed that the hydro- 
dynamic interactions become unimportant due to the effects of hydrodynamic screening [12]. To our 



^The data on page 5 and in Table 2 of this reference support our conjecture. 
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knowledge, screening has not been established in the theory of hard sphere suspensions. If it would 
occur in suspensions, the screened particle motion could be affected only by thermal fluctuations (truly 
Brownian motion!). Such Brownian hard spheres can be described by the short range interacting random 
walk model [13]. For finite concentrations, we expect the longer range hydrodynamic interactions to 
be very important. That this is indeed the case, is the central theme of our paper. In what follows, 
we provide the theoretical arguments in favor of hydrodynamic interparticlc interactions and screening 
which must be present in solutions at non-vanishing concentrations. By exploiting analogies between 
electrodynamics and fluid mechanics we shall demonstrate that hydrodynamic screening occurs in much 
the same way as screening of the magnetic field in superconductors. Therefore, mathematically, the de- 
scription of screening in suspensions is analogous to that for the Meissner effect in superconductors. This 
observation will allow us to account for a number of interesting properties of suspensions. For instance, 
the viscosity of hard sphere suspensions is known to diverge beyond some critical concentration ip* . This 
phenomenon has been observed experimentally and is well documented, e.g. see Refs. [14-18]. All these 
references arc concerned with changes in rlieological properties of suspensions occurring with changes in 
concentration (p. Using scaling-type arguments Brady, Ref.[19], and, independently, Bicerano et al, Ref. 
[20], had found that near = ip* the relative viscosity tj/ijq diverges as r]/r]Q = C(l — (p/(p*)~'^ with C 
being some constant. Furthermore, as it is shown by Bicerano et al, such analytical dependence of relative 
viscosity on concentration (p actually works extremely well for all concentrations. In view of (1.3), it is 
reasonable to expect vanishing of D for — > yj*. This phenomenon was indeed observed in Ref. [21]. . 

Theoretically, the result for relative viscosity was obtained as result of a combined nontrivial use of 
topological and combinatorial arguments. Such arguments can also be used, for instance, for description 
of the onset of turbulence in fluids or gases. As described in Ref. [22], such a regime in these substances 
is characterized by the sharp increase in the viscosity (just like in suspensions). According to Chorin, 
Ref. [22], Section 6.8, one can think about such an increase as analogous to processes which take place in 
superfluid ^He when one goes in temperatures from below to above A— transition, that is from the super- 
fluid to normal fluid state. Such a transition is believed to be associated with uninhibited proliferation 
of tangled vortices on any scale. In this work we demonstrate that Chorin's conjecture is indeed correct. 
This interpretation is possible only if both topological and combinatorial arguments are rigorously and 
carefully taken into account. Surprisingly, when this is done, the emerging description becomes isomorphic 
to that known for the Bose-Einstein condensation transition. Because of this, in addition to turbulence, 
in concluding section of this work we briefly discuss a number of apparently different physical systems 
whose behavior under certain conditions resembles that found in colloidal supensions. 

The rest of this paper is organized as follows. In section 2, we introduce notations, discuss experimental 
data with help of previously found generalized Stokes-Einstein relation [23] and make conjectures about 
how these results should be interpreted in the case if hydrodynamic screening does exist. Some familiarity 
with Ginzburg-Landau (G-L) theory of superconductivity is expected for proper understanding of this 
and the following sections. In Section 3 we study in detail how many particle diffusion processes should 
be affected if hydrodynamic interactions are taken into account. The major new results of this section 
are given in Section 3.3. where we rigorously demonstrate that account of hydrodynamic interactions 
causes modification of Pick's laws of diffusion in the same way as presence of electromagnetic field causes 
modification of the Schrodinger's equation for charged particles The gauge fields emerging in the modified 
Fick's equations are of zero curvature implying involvement of the Chern-Simons topological field theory. 
The following Section 4 considers in detail the implications of the results obtained in Section 3. The major 
new result of this section is given in Section 4.4. where we adopted the logic of the groimd breaking paper 
by London and London [24] in order to demonstrate the existence of hydrodynamic screening. Thus, the 
phenomenon of screening in suspensions is analogous to the Meissner effect in superconductors [25]. In 
Section 5 wc follow the logic of Ginzburg-Landau paper [26] elaborating the work by London brothers 
and develop similar G-L-type theory for suspensions. The major new result of this section is presented in 
Section 5.5. in which by using combinatorial and topological methods we reproduce the scaling results 
by Brady [19] and Bicerano et al, Ref. [20]. In Section 6 we place the obtained results in a much broader 
context. It is done with help of two key concepts: helicity and force-free fields. They had been in 
use for some time in areas such as magnetohydrodynamics, fluid, plasma and gas turbulence, classical 
mechanics written in hydrodynamic formalism but not in superconductivity or colloidal suspensions, etc. 
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In this section we mention as well other uses of these concepts in disciplines such as high temperature 
superconductivity, quantum chromo dynamics, string theory, non-Abelian fluids, etc. The paper also 
contains three appendices which are made sufficiently self contained. They are not only very helpful in 
providing details supporting the results of the main text but also of independentl interest. 



2 Stokes-Einstein Virial Expansions for a Broad Concentration 
Range 

2.1 General Results 

In 1976, Batchelor obtained the following general result for the cooperative diffusion coefficient [27]: 



K {(fi) ip / i9/i 



67r?7r 1 — \ dip 



(2.1) 



where K {p) is the sedimentation coefficient of the particles in suspension and is the chemical potential. 
Batchelor obtained for K {ip) the following result: 



K{<p) = 1-6.55(^ + 0(9?^) 



(2.2) 



so that (2.1) with thus obtained first order result for K [ip) can be used only for low concentrations. In 
Ref.[3] an attempt was made to extend Batchelor's results to higher concentrations. This was achieved 
in view of the fact that 



where 11 is the osmotic pressure. Use of this result in (2.1) produces: 



\ — p \ dip 



^^^^ QTrrjr (ydnj^j,' 



(2.3) 



(2.4) 



The Carnahan-Starling equation of state for hard spheres can be used to obtain the following result for 
compressibility 



an 

dn 

thus converting equation (2.4) into 



= kBT 



{i + 2^y + {p-4) ip3 



^(^)=f^w 

67r r/r 



(1 - p>r 



{l + 2pf + {ip-A) ip^ 



(2.5) 



(2.6) 



(1 - ip)' 

To be in accord with Batchelor's result (2.2) at low concentrations, the authors [3] suggested replacing of 
Eq.(2.2) by 

K ip) ^ {1 - p f-=' (2.7) 
which allows us to rewrite (2.6) in the following final form 



D/Do ^(l-p) 



6.55 



{i + 2py + {p>-4) ^3 



(2.8) 



(1 - pf 

convenient for comparison with experimental data. 

Such a comparison can be found in Fig. 12 of Ref [3] where this result is plotted against author's own 
experimental data for the cooperative diffusion coefficient. The experimental data within error margins 
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appears to agree extremely well with the theoretical curve obtainable from Eq.(2.8). However, it should 
be kept in mind that, in fact, originally Eq.(2.2) was determined only to first order in (and, therefore, 
only for the volume fractions less than about 0.05). Therefore, formally, Eq.(2.8) is in accord with Eq.(2.2) 
only for volume fractions of lesser than about 0.03. Therefore, it is clear from Fig 12 of [3] that to improve 
the agreement in the whole range of concentrations, a knowledge of a second order in is desirable in 
(2.2). This problem can be by passed as follows. 

Prom [3 ] the viscosity data from the same experiments were obtained so that the data can be fit to 
the following second order expansion: 

r?/r?o = 1 + 2.5(^ + 6.5V + O . (2.9) 

To obtain this result, the authors constrained the first order coefficient to 2.5 to comply with Einstein's 
result (1.2) for viscosity. If one considers these data without such a constraint, then one obtains, 

r?/% = 1 + 2.4^ + 7.1^2^0 ((^3^. (2.10) 

In the paper by Kholodenko and Douglas [23] the following result for the cooperative diffusion coeffi- 
cient was derived (the generalized Stokes-Einstein relation) 



D/Do 



5(0,0) 



So {0,0) 



-1/2 



(2.11) 



where S (0, 0) is the k = 0, zero angle static scattering form factor. The thermodynamic sum rule for the 
hard sphere gas produces the following result for this formfactor: 



5(0,0)T^/^ 



l + 4ip + 7ip'^ + 0(ip'^) . (2.12) 



5o(0,0). 

By combining Eq.s (2.9)-(2.12) the result for cooperative diffusion is obtained: 

D/Do = l + 1.6ip-3.9ip'^ + 0{ip^) . (2.13) 

For the sake of comparison with experiment, we made a numerical fit to the experimental data for higher 
concentrations obtained in [3] by a polynomial (up to a second order in <^) with the result*: 

D/Do = 1 + 1.5505<^ - 5.3663(^2 + O (if^) . (2.14) 

Comparison between Eq.s (2.13) and (2.14) shows that the theoretically obtained result, Eq.(2.13), is in 

good agreement with the experimental data, Eq.(2.14), within error margins. Alternatively, we can use 
the reciprocal of the empirical expression, Eq.(2.14), in (2.11) to obtain 

7?/% = 1 + 2.45(/j + 8.5(^2 + O ((^3^ , (2.15) 
which also compares well with the experimental data, Eq.(2.10). 

2.2 The generalized Stokes-Einstein relation and the role of hydrodynamic 
screening 

We return now to Eq.(2.11) for further discussion. Based on the results of introductory section, especially 
on Eq.s (1.1) and (1.3), we can formally write: 

D/Do = ^^ = l + a^ip + a2'P^ + 0{v^). (2.17) 
The actual values of , « = 1, 2, .. can be determined using Eq.(2.11) written in the following form 



*The correlation coefRcient obtained for this fit is 0.97. 
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(e/^o) 



^(0,0) 
5o(0,0) 



-1/2 



(2.18) 



in ho) 

where ^ is the correlation length, is the correlation length in the infinite dilution limit ~ R. To 
justify such a move we need to remind to our readers of some facts from the dynamical theory of linear 
response. To do so, we borrow some results from our previous work. [23] . 



Generally, both D and Do are measured by light scattering experiments. In these experiments the 
Fourier transform of the density-density correlator 

^(R, r) = {5n{r, t)Sn{r', t')) (2.19) 

is being measured. The formfactor, Eq.(2.19), is written with account of translational invariance, re- 
quiring the above correlator to be a function of relative distance |r — r'| = |R| only. Time homogeneity, 
makes it in addition to be a function of \t — t'\ =r only. In this expression, (...) represents an equilibrium 
thermal average while density fluctuations are given by 5n{v,t) = n{r,t) — {n). By definition, the Fourier 
transform of Eq.(2.19) is given by 



5(q,w) = dr dTS(R,T)e 



. Using this expression, we obtain the initial decay rate Fq ^ as follows [10,23]: 



J{<i-r-UT) 



2^ 



e--5(q,a;) 



With help of this result, the cooperative diffusion coefficient is obtained as 

In the limit of vanishingly low concentrations the self -diffusion coefficient is known to be [10] 

Do = iliml({r(t)-r(0)}^ 

b t—^00 t \ 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



where < ... > denotes the Gaussian- type average. Following Lovesey [28], it is convenient to rewrite this 
result as 

Dn = - I dt{v- v{t)) (2.24) 



in view of the fact that if 



c{t) - r(0) = j v{T)dT 



(2.25a) 



then, 




(r(i) - r(0))') = (I / v(r)dT I \ ^ 2 J dr J df (v(t) • v(f)) = 2 J drit-r) (v(r) • v(0)) , (2.25b) 



while, by definition, = | limt_>oo j ^(''(t) — r(0))' 

With these definitions in place and taking into account Eq.(2.18), wc would like now to discuss in 
more detail the relationship between D and Do- Using Eq.(2.29) of Ref. [23] we obtain 



D 



lim 

+0+,k: 



oo 



<j(0,t)-j(0,O> 
5(0,0) 



(2.26) 
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where the current j is given as j —dn{r, t)v(r, t), provided that the non-sHp boundary condition 



dr 



v(t) 



(2.27) 



is appUed. Here Vf{r,t) is the velocity of the fluid and r{t) is the position of the center of mass of the 
hard sphere with respect to the chosen frame of reference. Eq.(2.26) is in agreement with Eq.(2.24) in 
view of the fact that in the limit of zero concentration S{0, 0) = 1 so that < j(0, t) • j(0, t") >^ (v • v(t)) 
as we would like to demonstrate now. For this purpose, in view of Eq.(2.22) it is convenient to rewrite 
the result Eq.(2.26) is the equivalent form 



5(q,0)rW = j dt"ci- < j(q,t)j(-q,i") > -q |,^o+ 



(2.28) 



in accord with Eq.(2.15) of Ref.[23]. This relation is very convenient for theoretical analysis. For instance, 
it is straightforward to obtain D in the decoupling approximation as suggested by Ferrell [21]. It is given 

by 

q- < j(q, t)j(-q, t') > -q^q- < v(q', f)v(q - q', f') > •q(5n(q - q', t)Sn{q' , t')). (2.29) 

In Section 5.4. we provide proof that the above decoupling is in fact exact. This provides an explanation 
why it is working so well in real experiments. 

In the meantime, wc shall consider this decoupling as an approximation. Once such an approximation 
is made, the problem of calculation of D is reduced to the evaluation of correlators defined in Eq. (2.29) . For 
the velocity-velocity correlator, the following expression was obtained before (e.g. see Ref.[23], Eq.(2.18)): 



< Vii,{t)vji,>{t') >= (27r)3(5(k + k'){(5y - - 



fc2 r]k'^ 



= 2kBTnij{\s.)S{t-t') 



(2.30) 



with i,j = l — 3. This expression defines the Oseen tensor Hij{k) to be discussed in detail in the next 
section. The presence of the delta function 5{t — t') in Eq.(2.30) makes it possible to look only at the 
equal time density-density correlator in the decomposition of the j — j correlator given by Eq.(2.29). Such 
a correlator also was discussed in Ref.[23] where it is shown to be 



((5n(k - k', t)5n{\s!, t)) = ksT <n> 



d_ 



<n> 



(2.31a) 



Actually, it is both time and fc— independent since, as is well known, it is the thermodynamic sum rule. 
That is 

" d 

5(0,0) = ksT <n> 



dIL 



<n> 



(2.31b) 



It is convenient to rewrite this result as follows 



knT < n > 



an 



<n> 



dRS(R,0) =5(0,0), 



implying that 



SCRfi) =kBT <n> 



an 



< n > 



(5(R). 



(2.31c) 



(2.32) 



In view of Eq.(2.12), we notice that in the limit of vanishing concentrations 5(0, 0) = l.In such an extreme 
case the decoupling made in Eq.(2.29) superimposed with the definition, Eq.(2.22), and the fact that 



d 1 d d 



dqi dqj 
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produces the anticipated result, Eq.(2.24), as required. Next, following Ferrell [21], we regularize the 

CO 

delta function in Eq.(2.32). Using an identity 1 = / dxxe~'^, the regularized expression for 5(R,0) is 



obtained as follows 



5(r,0) 



_kBT 



<n> 



d_ 

m 



<n> 



(2.33) 



where r = |R| and the parameter ^ is proportional to the static correlation length^. To use this expression 
for calculations of D, by employing Eq.(2.26) we have to transform the hydrodynamic correlator, Eq.(2.30), 
into coordinate form as well. Such a form is given in Eq.(2.33) of Ref.[23] as 



<v{r,t)-v{r',t') >= 



keT 1 



TTT] |r — r 



-Mt-t'). 



(2.34) 



This expression is written with total disregard of possible effects of the hydrodynamic screening, though. 
The combined use of Eq.s (2.33) and (2.34) in Eq.(2.26) produces the anticipated result 



D 



37r?7^ 



(2.35) 



in accord with that obtained by Ferrell, Rcf.[21], Eq.(ll), provided that we redefine (still arbitrary) the 
parameter ^ as 2^. The result (2.35) also coincides with Eq.(1.3) if we identify ^ with R*. Furthermore, 
by looking at Eq.(2.18) wc realize that in the infinite dilution limit we have to replace ^ by and, 
accordingly, 77 by ?7g. Such an identification leads to the generalized Stokes-Einstein relation in the form 
given by Eq.(2.18) implying that 

5(0,0) T'/^ 



5o(O,0). 

Since we have noticed before that (0, 0) = 1 this result can be rewritten as 



(2.36a) 



(2.36b) 



Suppose now that hydrodynamic interactions are screened to some extent. In such a case the result 
Eq.(2.34) should be modified accordingly. Thus, we obtain 



< v(r,f) • v(r',t') >= 



kBT 



exp(-— ) 



(2.37) 



where we have introduced the hydrodynamic correlation length ^jj. If, as we shall demonstrate below, the 

analogy between hydrodynamic and superconductivity makes sense under some conditions then, using 
this assumed analogy we introduce the Ginzburg parameter kq for this problem via known relation [25]: 

= i^gI (2.38) 
Using Eq.s (2.33), (2.37) and (2.38) in (2.26), the result for D, Eq.(2.35), acquires the following form: 

1 



67rSj?^ 



+ — )" 



KG 



(2.39) 



Since, the adjustable parameter S is introduced in Eq.(2.39) quite arbitrarily, we can, following Ferrell, 
Ref. [21], take full advantage of this fact now. To do so, we notice that from the point of view of the 
observer, the relation (2.36) holds irrespective of the absence or presence of hydrodynamic screening. 
Because of this, we write 

S(l + — )=^ (2.40) 
Kg 



^For more details, see our work, Ref. [23]. 
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so that the Eq.(2.36) used in the generaUzed Stokes-Einstein relation remains unchanged. By combining 
Eq.s(2.36b), (2.38) and (2.40) we obtain: 

KG = ^ , ^ = # (1 + — ) (2.41a) 

or, equivalently, 

^= ,^ =(! + —)■ (2.41b) 

S (0,0)^ ^ ^ 

In this equation the parameter E is still undefined. We can define this parameter now based on physical 
arguments. In particular, let us set S = S (0, 0) ^q. Then, we end up with the equation 

1 + — = , ^ (2.42) 

KG (0,0) 

leading to 

KG = 1 ^ ^ - (2.43) 
To reveal the physical meaning of this equation we use Eq.s (2.36b), (2.38) and (2.43) in order to obtain 

= ^i^. ,2.44, 



From Eq.(2.12) we notice that by considering the infinite dilution limit 0, we obtain: oo, 

implying absence of hydrodynamic screening. Consider the opposite case: ^ oo (that is, in practice, 
being large). Looking at Eq.(2.12) we notice that in this case ^ ^ indicating the complete screening, 
as expected. Using Eq.(2.42), these results allow us to rewrite the generalized Stokes -Einstein relation, 
Eq.(2.18), in the equivalent form emphasizing the role of hydrodynamic screening. Thus, we obtain, 

£>/£>o = ^^(l + — ). (2.45) 
(»?/%) kg' 



3 Diffusion processes in the presence of hydrodynamic interac- 
tions 

3.1 Some facts from the diffusion theory 

If n is the local density, then the flux j= nv obeys Pick's first law: 

j = -DVn, (3.1) 

where D is the (in general, cooperative) difi^usion coeSicient, and v is the local velocity. Upon substitution 
of this expression into the continuity equation 

f + V.j = (3.2) 



we obtain the diffusion equation commonly known as Pick's second law 

fin 

= DV^n. (3.3) 



dt 
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In the presence of some external forces, i.e. 



F = -VC/, (3.4) 

the diffusion laws must be modified. This is achieved by assuming the existence of some kind of friction, 
i.e. by assuming that there exists a relation 

7V = F (3.5) 

between the local velocity v and force F with the coefficient of proportionality 7 being, for instance (in 
the case of hard spheres), of the type given in Eq.(l.l). With such an assumption, the diffusion current, 
Eq.(3.1), is modified now as follows 

j = -DVn - -VU. (3.6) 

7 

Such a definition makes sense. Indeed, in the case of equilibrium, when the concentration Ugq obeys the 
Boltzmann's law ^ 

Ueq = noexp(-— -), (3.7) 

KB-L 

vanishing of the current in Eq.(3.6) is assured by substitution of Eq.(3.7) into Eq.(3.6) thus leading to 

the already cited Einstein result. Eq.(l.l), for Dq. As in the case of Eq.(1.3), we shall assume that for 
finite concentrations one can still use the Einstein-like result for the diffusion coefficient. With such an 
assumption, the current j in Eq.(3.6) acquires the following form [12]: 

n Th 

j = V{kBT\nn + U) = V/x, (3.8) 

7 7 

where the last equality defines the nonequilibrium chemical potential fx, e.g. like that given in Eq.(2.1). 
Alternatively, the modified flux velocity v/ is given now by —^'Vfj, so that the continuity Eq.(3.2) reads 
as 

^ + V.(nv^) = 0. (3.9) 

Exactly the same equation can be written for the probability density ^ if we formally replace n by \1/ 
in the above equation [10]. Such an interpretation of diffusion is convenient since it allows one to talk 
about diffusion in terms of the trajectories of Brownian motion of individual particles whose positions 
Xn{t),n = 1,2,... are considered to be as random variables. Then, the probability ^ describes such 
collective Brownian motion process described by the following Schrodinger-like equation 

n 

in which the velocity v/„ is given by 

V/n = -^inm^(fcBTln* + C/). (3.11) 

Thus, we obtain our final result 

= ,3.1.) 

m.n 

adaptable for hydrodynamic extension. For this purpose, we need to remind our readers of some basic 
facts from hydrodynamics 
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3.2 Hydrodynamic fluctuations and Oseen tensor 

The analog of Newton's equation for fluids is the Navier-Stokes equation. It is given by [29] 

^v + (v V)v = -^VP + rvV (3.13) 

where P is the hydrodynamic pressure, F = ilo/Po the kinematic viscosity and Pq is the density of the 
the pure solvent. At low Reynold's numbers, the convective term (v-'V)v can be neglected [29], p. 63. We 
shall also assume that the fluid is incompressible, i.e 

div V = 0. (3.14) 

Under such conditions the Fourier transformed Navier-Stokes equation can be written as 

^Vk = -rkV--Pk. (3.15) 

ot p 

Let us add a fluctuating source term fk to the right hand side of Eq.(3.15). Then, using the incompress- 
ibility condition, Eq.(3.14), we obtain: 

a = -ip ^2 • (3-16) 
Introducing the transverse tensor Tjj(k) = 5ij — and decomposing a random force as 

fl{t) = Y.T,i{\^)f^^{t) (3.17) 

3 

eventually replaces the Navier-Stokes equation by the Langevin-type equation for the transverse velocity 
fluctuations: 

^Vk + rk2vk = fJ(i)- (3.18) 

As is usually done for such type of equations, we shall assume that the random fluctuating forces are 
Gaussianly distributed. This assumption is equivalent to the statement that 



[flmfu'it )) = T,j{k){2nfS{k + k')D6{t - t') (3.19) 
with parameter D to be determined. A formal solution of the Langevin-type Eq.(3.18) is given by 



t 



Mt) = Vk(0)e-" + J dt'fl it') e-f (*-*') (3.20) 



with r = k^r. Introducing Vk(t) — Vk(0)e~'"* = Vk(t), we obtain 

, t t' , 

{vi^m^'it')) = / j dt'i^ it') e-^(*-*') I dt"tl{t")e-^^''-'"^ \ . (3.21) 
^0 ' 

To calculate this correlator, and to determine the parameter 2D, we consider the equal time correlator 
first. In such a case the equipartition theorem produces the following result: 

{Vi^,{t)v^y,,{t)) = r,,(k)(27r)3^(k + k')^. (3.22) 
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Taking into account Eq.s (3.19) and (3.22) we obtain for the velocity-velocity correlator, Eq.(3.21), the 
following result 

{vMt)vjk'it')} = Ti,(k)(27r)35(k + k')^Jexp(-f |t-t'|) 



pT 2 

r 

'2 



= 2fcBm,(k)-exp(-f|t-i'|). (3.23) 



In the limit f — > 00 the combination ^ cxp(— f |t — t'|) can be replaced by S{t — t'). In this limit the 
obtained expression coincides with already cited Eq.(2.30). Furthermore, the constant D can be chosen 
as To prove the correctness of these assumptions, we take a Fourier transform (in time variable) in 
order to obtain 

{vM %(-k)(-a;)) =^^^T,,(k). (3.24) 

This result coincides with Eq.(89.17) of Ref.[25] as required. Here the sign = means "up to a delta function 

prefactor" . Incidentally, these prefactors were preserved in another volume of Landau and Lifshitz, e.g. 
see Ref. [30], Eq.(122.12). Since in the limit a; ^ we reobtain (upon inverse Fourier transform in time) 

Eq.(2.30), this fact provides the needed justification for replacement of the factor ^ exp(— F |t — t'|) by 

5{t-t'). 

In polymer physics, Ref. [10], typically only this w ^ limit is considered, which is equivalent to 
considering physical processes at time scales much larger than the characteristic time scale = pB^/nrjQ 
mentioned in the Introduction. Although this fact could cause some inconsistencies (e.g. sec discussion 
below), we shall follow the traditional pathway by considering mainly this limit causing us to drop 
altogether time-dependence in Eq.(3.15) thus bringing it to the form considered in the book by Doi and 
Edwards, Ref. [10], Eq. (3. III. 2). Following these authors, this approximation allows us to specify a 
random force f(r) as 

f (r) = J2 FnS{r - R„) (3.25) 

n 

implying that particle (hard sphere) locations are at the points R„ so that the fluctuating component of 
fluid velocity v(r) at r is given by 

v(r) =^H(r-R„)-F„ (3.26) 

n 

with the Oseen tensor H^j (r) in the coordinate representation given by 

7* ' 

In this expression fi = yj-^^ view of Eq. (2.27), we can rewrite Eq.(3.26) in the following suggestive form 

v(R„)= H(R„-R„)-F„, (3.28) 

for velocity v(R„) of the particle located at R„. 

3.3 Fick's laws in the presence of hydrodynamic interactions. Emergence of 
gauge fields 

By comparing Eq.s(3.12) and (3.28) we could write the Fick's first law explicitly should the Oseen tensor 
be also defined for m = n. But it is not defined in this case. As in electrostatics, self-interactions must 
be excluded from consideration. In view of the results of Section 2, the situation can be repaired if we 
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assume that at some concentrations the hydrodynamic interactions are totally screened. In such a case 
only the usual Brownian motion of individual particles is expected to survive. With these remarks, Pick's 
first law for such hydro dynamically interacting suspensions of spheres can be written now as follows 

v/(R„) = - Vh(R„ - R„) • -|-(fcBTln* + U), (3.29) 

where the redefined Osccn tensor Hy (R) has the diagonal part Hii(R) = ^ in accord with Eq.(l.l). The 
potential U comes from short-range non- hydrodynamic interactions between particles, which are always 
present. Using this result and Eq.(3.10), we finally arrive at the Pick's second law 



n.m 



in accord with Eq.(3.110) of Ref.[10]. Since this equation contains both diagonal and nondiagonal terms 

the question arises about its mathematical meaning. That is, we should inquire: under what conditions 
does the solution to this equation exist? The solution will exist if and only if the above equation can be 
brought to the diagonal form. To do so, as it is usually done in mathematics, we have to find generalized 
coordinates in which the above equation will acquire the diagonal form. Although the attempts to do so 
were made by several authors, most notably, by Kirkwood, e.g. see Ref.[10], chr-3 and references therein, 
in this work we would like to extend their results to account for effects of gauge invariance. 

We begin with the following auxiliary problem: since =div- V = V • V, wc are interested in finding 
how this result changes if we transform it from the flat Euclidean space to the space described in terms 
of generalized coordinates. This task is easy if we take into account that in the Euclidean space 

with /i'-' being a diagonal matrix with unit entries. We notice that the above expression is a scalar and, 
hence, it is covariant. This means, that we can replace the usual derivatives by covariant derivatives, the 
metric tensor h^^ by the metric tensor g'^^ in the curved space so that in this, the most general case, we 
obtain 

D,g'W,f{^) = A5^._^/(x) +/ir^, A/(x), (3.32) 

where summation over repeated indices is assruned, as usual. The covariant derivative Dj is defined for 
a scalar f as Dif = and for contravariant vector X' as 

DiX^ = ^ + T)^X^ (3.33) 

with Christoffel symbol F^-^ defined in a usual way of Riemannian geometry. A precise definition of this 
symbol is going to be given below. Since T\^. = -^In^ , we can rewrite Eq.(3.32) in the following 
alternative final form 

VV = A.'^I^./(x) = -^^[,^.V^^/] (3.34) 

so that in Eq.(3.3) the operator is replaced now by that given by Eq.(3.34). To make this presentation 
complete, we have to include the relation 

In the simplest case, when we are dealing with 3 dimensional vectors, so that r = r((j'i, g2, 93), sometimes 
it is convenient to introduce vectors 

= (3-36) 
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and the metric tensor 

Qij = e, • Bj (3.37) 

with "•" being the usual Euclidean scalar product sign. Definitions Eq.(3.35) and (3.37) are obviously 
equivalent in the present case. Because of this, it is clear that upon transformation to the curvilinear 
coordinates the Riemann curvature tensor written in terms of gij is still zero since it is obviously zero for 
the hki- The curvature tensor will be introduced and discussed below. Before doing so, using the example 
we have just described, we need to rewrite Eq.(3.30) in terms of generalized coordinates. In the present 
case, we must have 3A'' generalized coordinates and the tensor hki is not a unit tensor anymore. Our 
arguments will not change if we replace Eq.(3.30) by that in which the potential U = 0. Furthermore, 
we shall adsorb the factor fc^T into the tensor H and this redefined tensor we shall use instead of hki- 
Evidently, the final result for the Laplacian, Eq.(3.34), will remain unchanged. The question arises: if in 
the first example the Riemannian curvature tensor remains flat after a coordinate transformation (since 
the tensor hki is the tensor describing the flat Euclidean space), what can be said about the Riemann 
tensor in the present case? To answer this question consider once again Eq.(3.35), this time with the 
tensor H„i„ instead of h^j. For the sake of argument, let us ignore for a moment the fact that each of 
generalized coordinates is 3-dimensional. Then, we obtain, 

dR''- dR^ ,^ ^„ , 

where we introduced a set of new generalized coordinates {Q} = {gi,...,giv} so that Ri = Ri{{Q}). We 
shall use Greek indices for new coordinates and Latin for old. In the case of 3 dimensions the above result 
becomes 

^?°/3 = ^ • Hki • ^ (3.38b) 

with "•" being the Euclidean scalar product sign as before. The indices k, 1, a and /3 now have 3 
components each. We are interested in generalized coordinates which make the metric tensor go,p diagonal. 
By analogy with Eq.(3.36), we introduce now a scalar product 

< R • R >= R"" • Hki • R' (3.39) 

so that instead of the vectors e, we obtain now 

and, accordingly, 

gap =< e„ • > . (3.41) 
The Christoffel symbol can be defined now as 

To find the needed generalized coordinates, we impose an additional constraint 

9e„ dep 



(3.43) 



compatible with the symmetry of the tensor Hki. By combining Eq.s(3.42) and (3.43) we obtain, 

r>^=r^„e^ (3.44) 

implying that — F^^. That is, the imposition of the constraint, Eq.(3.43), is equivalent to requiring 
that our new generalized spacc^ is Riemannian (that is, without torsion). In such a space we would like 
to consider the following combination 
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Again, using Eq.(3.42) we obtain 

a^(rl.e,) = (^r^,) + (^e,) Tl,. (3.46) 



Analogously, we obtain 

^(^e.) = (^a) + (^e,) r^„. (3.47) 
Finally, we use Eq.(3.42) in Eq.(3.46) and (3.47) in order to obtain the following result for Rap 

= (3-48) 

The second line defines the Riemann curvature tensor. In the most general case it is given by R^g/^- By 
combining Eq.s(3.40), (3.43), (3.45) and (3.48) we conclude that 

- ^'^^ (3.49) 



implying that the Riemann tensor is zero so that the connection is flat. For such a case we can 
replace the covariant derivative Di by Vj +Aj [31]. The vector field Ai is defined as follows. Introduce a 
1-form A via A = Aidx'^, Ai = AfT", where in the non -Abelian case is one of infinitesimal generators 
of some Lie group G obeying the commutation relations [T",T^] — if^^^T^ of the associated with it 
Lie algebra. In addition, tr[T«r'^] = ^6"'^. The Chern-Simons (C-S) functional CS{A) producing upon 
minimization the needed flat connections is given by[31,32] 

CS{A) = ^ [ tr{A A dA + |a A A A A) 
47r J 3 

M 

= -^j s^'^triA^idjA, - dkAj) + lA,{A„Ak]) (3.50) 



M 



with k being some integer. Minimization of this functional produces an equation for the flat connections. 
Indeed, we have 

y (^^(A + B) = j tr{B ArfA + AAB + 2BAAAA) 

M 

= 2 j tr{BA{dA + AAA), (3.51) 

M 

where we took into account that 

/BR r f) A 

tr{Aidx' A g^dx^ A dx'') = / tr{Bkdx^ A Q^dx' A dx'). 

M M 

From here, by requiring 

— CS(A + B) = (3.52) 

on 

we obtain our final result: 

OA- dA- 

dA + AAA={—^-—^ + [Ai, Aj])dx' A dx^ = F{A)dx' A dx^ = 0. (3.53) 

(-/Jbj (J Jb ^ 

In the last equality we have taken into account that both in the C-S and Yang-Mills theory F{A) is 
the ciu'vature associated with connection A. Vanishing of curvature produces Eq.(3.53) for the field A. 
Irrespective to the explicit form of the field A, we have just demonstrated that, at least in the case when 
the potential U in Eq.(3.30) is zero, this equation can be brought into diagonal form provided that the 
operator V, is replaced by Vj +Ai with the field Ai to be specified below, in the next section. 
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4 An interplay between topology and randomness: connections 
with the vortex model of superfluid ^He 



4.1 General comments 



The C-S functional, Eq.(3.50), whose minimization produces Eq.(3.53) for the field A was introduced 
into physics by Witten [32] and was discussed in the context of polymer physics in our previous works 
summarized in Ref.[33]. Since polymer physics of fully flexible polymer chains involves diffusion- type 
equations [10], the connections between polymer and colloidal physics are apparent. For this reason, we 
follow Ref. [32] in our exposition and use it as general source of information. 

Specifically, as explained by Witten [32], theories based on the C-S functional are known as topological 
field theories. The averages in these theories produce all kinds of topological invariants (depending upon 
the generators T" in the non Abclian case) which arc obscrvablcs for such theories. In the present case 
the question arises: should we use the non Abelian version of the C-S field theory or is it sufficient to 
use only its Abelian version, to be defined shortly? Since both versions of C-S theory were discussed in 
the context of polymer physics in Ref. [33], wc would like to argue that, for the purposes of this work, the 
Abelian version of the C-S theory is sufficient. We shall provide the proof of this fact in this section. 

The action functional for the abelian C-S field theory is given by® 

S^_s[A] = d^xe^^PA^d^A, (4.1) 

M 

With such defined functional one calculates the (topological) averages with help of the C-S probability 
measure 

< • • • >c-s= N j D[A] exp{i5^_s[A]} • • • . (4.2) 
The random objects which are subject to averaging are the Abelian Wilson loops W{C) defined by 

W{C) = exp{ie j> dr ■ A}, (4.3) 

c 

where C is some closed contour in 3 dimensional space (normally, without self-intersections) , and e is some 
constant ("charge") whose exact value is of no interest to us at this moment. The averages of products 
of Wilson's loops (perhaps, forming a link L) 

n 

W{L) = l[W{Ci) (4.4) 

'i=i 

are the main objects of study in such a topological field theory. Substitution of W{L) into Eq.(4.2) 
produces the following result [32] 

< W{L) >c-s= exp{i (y ) E (^i^Mh J)} (4-5) 
with the (Gauss) linking number lk(i,j) defined as 

T T 

1 iiu„ (^i-r,-) _ 1 f i (r(si) - r(s,)) 



Hhj) = ^ff ^'^^i ^ '^''i] ■ I _ ,3 = ^ I I dsidsj[v{si) X v(sj)] • -^-j—^ — Z7~^- (^"^^ 



C, Cj 



|r(si) - r(sj)|^ 



^E.g. see Eq.(4.12) ofRef.[33]. 
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Here Tj and Tj are respectively the contour lengths of contours Cj and Cj and v(s) = j^r{s). With 
the Gauss linking number defined in such a way, in view of Eq.(4.5), it should be clear that we must to 
consider as well self-linking numbers lk{i, i). Such a technicality requires us to think about the so called 
framing operation discussed in some detail in both Ref.s [32] and [33]. We shall ignore this technicality 
until Section 6 for reasons which will become apparent. 



4.2 An interplay between the topology and randomness in hydrodynamics 

Following Tanaka and Ferrari, Refs [34, 35], we rewrite the Gauss linking number in a more physically 
suggestive form. For this purpose, we introduce the "magnetic" field B(r) via 

B(r) = — / dr, X (4.7) 
4^7 |r-r(s,-)| 

allowing us to rewrite the linking number lk(i,j) as 

lk{i,j) = j>dVi-B{vi). (4.8) 

Ci 

Eq.(4.7) for the field B(r) is known in magnetostatics as Biot-Savart law, e.g. see Ref.[36], Eq.(30.14). 
Because of this, we recognize that 

V • B =0 (4.9) 

and 

VxB=j, (4.10) 

where j(r) = § ds\'{s)5{Y — r(,s)).To connect these results with hydrodynamics, we introduce the vector 

c 

potential A in such a way that 

VxA = B. (4.11) 
Using this result in Eq.(4.10) we obtain as well 

V'A = -j (4.12) 

in view of the fact that V • A =0. In hydrodynamics we can represent the local fluid velocity following 
Ref.[37], page 86, as 

V = V X A (4.13) 

and define the vorticity u} as 

(I; = V X V. (4.14) 
By analagy with Eq.s(4.10) and (4.12) we now obtain 

V^A = -w. (4.15) 

Hence, to apply previous results to hydrodynamics the following identification should be made: 

u^Uv^'B. (4.16) 

The kinetic energy f of a fluid in a volume M is given by 

_P 

M 



£ =^ / v^d^r. (4.17) 



We would like now to explain how this energy is related to the above deflned linking numbers. For this 
purpose, we introduce the following auxiliary functional: 



(4.18) 
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Use of the theorem by Stokes produces 

jr[A], = j^dvi- A(ri) = jj dSi-i-^ xA) = JJdSi-w=JJdSi^ dsjv{sj)S{r.-r{sj)). (4.19) 

C, 5, Si Si Cj 

At the same time, for the hnking number, Eq.(4.8), an analogous procedure leads to the following chain 
of equalities 

lk{i,j) = jdvi- B(rj) = jdvi- v{Yj) = jj dSi-{V xv)= jj dSi-uf. (4.20a) 

Ci Ci Si Si 

Since the same vector potential was used in both Eq.s(4.11) and (4.13) we notice that Eq.s (4.12) and 
(4.15) also imply that 

u) =ej, (4.20b) 
where e is some constant to be determined. Because of this, we obtain 

eJ^[A]i = lk{i,j) = eJJ dS, dsjv{sj)S{r,-r{sj)). (4.21) 

S, Cj 

Since the obtained equivalence is of central importance for the entire work, we would like to discuss a few 
additional details of immediate relevance. In particular, from Eq. (4.20b), which we shall call from now 
on, the London equation (e.g. sec the Subsection 4.4 below), it should be clear that the as yet imknown 
constant e must have dimensionality of inverse length L~^. This fact should be taken into account when 
we consider the following dimensionless'' functional 

^t^^ = 2^ J '^''■^'^ ^ + 7 E / '^^^ ■ (4-22) 

M ^ Cj 

and the path integral associated with it, i.e. 

N j D[A]5{V ■ A) exp{-W[A]} =< W{L) >t (4.23a) 
to be compared with Eq.s(4.2) and (4.5). Here the thermal average < • • • >t is defined by 

<...>T=N j D[A]5{V ■ A) exp{-^^ j d^r{V x A)'} • • • . (4.23b) 



M 



Calculation of this Gaussian path integral is complicated by the presence of a delta constraint (Coulomb 
gauge) in the path integral measure. Fortunately, this path integral can be found in the paper by 
Brereton and Shah [38] . Without providing the details, these authors presented the following final result 
in notations adapted to this work: 



2 * * 



< W{L) >T= exp{-^ (y) ^. .^^ J J dSidsjr{s,) ■ H[r(s,) - r(s,)] • f (s,)}. (4.24) 



The Oseen tensor H(R) in this expression was previously defined in Eq.(3.27) and the prime on the 
summation sign means that the diagonal part of this tensor should be excluded. Even though calculations 

leading to this result are not given in Rcf. [38] , they can be easily understood ficld-thcoretically. For this 
purpose, wc have to regularize the delta function constraint in the path integral measure in Eq.(4.23) very 



'^In view of the fact that the dimensionahty of e is fixed we have introduced a factor / which makes the functional W[A] 
dimensionless. This factor will be determined shortly below. 



18 



much the same way as Ferrell, Ref.[21], did it in the case of hydrodynamics as we discussed in Section 2. 
Specifically, we write 

N J D[a{r)] exp(-^ J Sva^{v)) J D[A]d{V • A - a{r)) exp{-^ J d^r{V x A)^} • • • 

M 

= nJ D[A] exp{- 2^ I d^V X A)^ - 1 1 rfV(V • A)^} • • • 

M 

= nJ D[A] exp{-^ J dhA^[-5^^V^ - (1 - l)d^d^]A^} (4.25) 



with some adjustable regularizing parameter ^.Also, for the quadratic form (in A) in the exponent of the 
last expression we obtain 

/ d^rA^[-5^,V^ - (1 - l)d^d,]A, = I rf3kA^(k)[,5^,k2 - (1 - l)k^K\A, (4.26) 

The inverse of the matrix [iS^j^k^ ~ (1 — ^)kfikv\ is easy to find following Ramond [39]. Indeed, we write 

[5^X - (1 - l)k^K][X{\<i)5.p + y(k)fc.fc^] = 5^p. (4.27) 

From here the unknown functions ^(k) and l'(k) can be determined so the inverse matrix is given 
explicitly by 

[X{\.)5., + y(k)fc.M = ^[5.p - (1 - l)^]. (4.28) 

In the limit ^ — > wc recover the Oseen tensor (up to a constant l/r;) in the A;-space representation in 
accord with Ref. [10]. These results explain why in the average, Eq.(4.24), there are no diagonal terms. 
Now we are ready to determine the constant e introduced in Eq.(4.22). 

4.3 Repcirametrization invciriance and vortex-vortex interactions 

The important result for < W{L) >t contains random velocities f (s) and thus, seemingly, additional 
averaging is required. The task now lies in finding the explicit form of this averaging. To do so, several 
steps are required. To begin, we notice that in the absence of hydrodynamic interactions Eq.(3.30) 
acquires the following form 

with diS'usion coefficient Dq defined in Eq.(l.l). If Eq.(4.29) we treat \1/ as Green's function (e.g. see 
Appendix A for details), then it can be formally represented in the path integral form as 

t 

^{t; Ri,... Rn) = I D[{K{t)}] exp(-^ J E^_^ [f (r,)]' dr,). (4.30) 
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In this expression we have suppressed the explicit R-dependence of the path integral to avoid excessive 
notation. Hydrodynamic interactions can now be accounted for as follows 



^ = nJ D[A] exp{-^ J d\A^[-d^.V^ - (1 - hd^dM 

M 

t t 

X J Z)[{r(r)}]exp(-^|(5]J_Jf(r,)]^dr,)exp{z^ J J^(r,)] • A[r(r,)]dr,} 

' 

N * 

= < n / DMrj)}] e^P(-i^ / dr^) exp{z^ j r(r,) • A[r(r,)]dT,} >^. (4.31) 

^'=^ 

Perturbative calculation of path integrals of the type 

t t 

I[A; t]=J D[{r{T^)}] exp(-^ J [^(r,)]' drj exp{i^ J [r(r,)] • A[r(r,)]dr,} (4.32) 



was considered by Feynman long ago, Ref. [40] . From this paper it follows that the most obvious way to 
do such a calculation is to write the usual Schrodinger-like equation 

(J-^ - Do{P - ieAf^ G{t, r; t'r') = 0, r ^ r' (4.33) 

and to take into account that (P — ieA)^ — F^ — ieA ■ P—ieP ■ A — e^A^ ~ — ieA ■ P+0(A^) (since 
P • A =0). This result is useful to compare with Eq.(3.32) in order to recognize that the field A is indeed 
a connection. 

To use these results, we would like to rewrite Eq.(3.30) in the alternative form which (for [/ = 0) is 
given by 

In arriving at this equation we took into account Eq.(3.14). Consider such an equation for n = 2. In this 
case, we rewrite Eq.(4.34) in the style of quantum mechanics, i.e. 

^ - ifi - ^2 - * = (4.35) 

in which, as in quantum mechanics, we shall treat V\2 as a perturbation. The best way of dealing with 
such problems is to use the method of Green's functions. For our reader's convenience we present some 
facts about this method in Appendix A. Eq.(A.lO) of this Appendix provides an equation for the effective 
potential V. A similar type of equation was obtained in the book by Doi and Edwards, Ref. [10], in Section 
5.7.3., who used methods of the ofFcctivo medium theory. Using this theory they were able to prove 
the existence of screening for the case of polymer solutions. We shall reach an analogous conclusion 
about screening in colloidal suspension using different arguments to be discussed in the next subsection. 
In the meantime, we would like to provide arguments justifying our previously made approximation: 
(P — ieA)"^ ~ P^ — ieA ■ P+0(A^). Using results of Appendix A, we introduce the one-particle Green's 
function Go as a solution to the equation 

[jt ~ ^^^2) <^o(R"*;R"*') = ^(R-- - *')• (4-36) 

Having in mind the determination of previously introduced factor / (in Eq.(4.22)), it is convenient to 
rescale the variables in this equation to convert it into dimensionless form. Evidently, the most convenient 
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choice is t = t/{Dq/Rq) and R = i?oR. with i?o is the hard sphere radius introduced in Eq.(l.l) and r 
and R being dimensionless time and space coordinates. Below, we shall avoid the use of tildas for R and 
shall still write t instead of r.The original symbols can be restored whenever they are required. Having 
this in mind, next we consider the two-particle Green's function Go -In the absence of interactions, it is 

just the product of two Green's functions of the type given by Eq.(4.36). As a result, the Dyson-type 
equation for the full Green's function for Eq.(4.34) (n = 2 case) is given by 

G(Ri, R2, Ri , R2 , t") = Go(Ri, R'l, <:')Go(R25 R2> i') 

-j- j (jo(Rl) ^! Rl) ^l)^o(R2) t't R2? ^l)^(-^l5 ■^2)^(-^l^ -^2) ^1) -^1 5 -^2 ? ^")c^Ri^^R2^^1 

(4.37) 

in which the potential ^(R'jjRj) = ksTiiijl'Ki — R2) aij"' before, summation over repeated 
indices is assumed. Using results of Appendix A and Eq.(4.37) it is possible to write now the equation 
for the effective potential. In view of the results to be discussed in the next subsection, this is actually 
unnecessary. Hence, we proceed with other tasks at this point. Specifically, taking into account Eq.(3.27) 
in which the explicit form of the Oseen tensor is given, we conclude that the nondiagonal part of this 
tensor can be discarded in the Dyson Eq.(4.37). This is so because of the following obvious identity: 
[(ri — r2) • ri] [(ri — r2) • r2] + [(r2 — ri) • ri] [(r2 — ri) • r2] =0 associated with the scalar products of 
unit vectors in Eq.(3.27). Evidently, it is always possible to select a coordinate system centered, say, 
at ri Alternatively, this result can be easily proven in A;-space taking into account the incompressibility 
constraint. Furthermore, these observations cause us to write the potential y(Ri,R2) in the following 
dimensionful form^ 

Using dimensional analysis of Eq.(4.36), this result can be easily rewritten also in dimensionless form. 
Explicitly, it is given by 

^C''-^' - 4.X^Do|R.-R.| as; ^ C^'"") 

in which the scalar product can be of any sign. This fact is of importance because of the following. 

Using Eq.(4.31) and proceeding with calculations of the path integral following Feynman's prescrip- 
tions [40], we obtain exactly the same equation as that given by Eq.(4.37). This observation allows us to 

determine the constants e and / explicitly. In view of the results just obtained, the constant e can be 
determined only with accuracy up to a sign. Taking this into account, the value of e is determined as 

e = ± — 1/ — —, while the constant f is given by Dr. in view of the fact that the field A in Eq.(4.22) has 
i?o V 47r?7 

dimensionality LF' /t , i.e. that of the diffusion coefficient, while the dimensionality of e is fixed by the 
Eq. (4.20b), so that the combination erfsr(s) is dimensionless. 

Using these results and Eq.(4.38), we can rewrite < W(L) >t defined by Eq.(4.24) in the following 
manifestly dimensionless physically suggestive form 

ksT Si Sj / / \dr{Ti) ■ dr{Tj)\ . 



< W(L) >T= exp(--^ii^ > ZLZ^ \JJ1\ (4.39) 

^ DoStttj ^i,3=i R„ R„ J J |r(T,)-r(T,)| ^ ^ ^ 

where we have introduced the dimensionless Ising spin-like variables Sj playing the role of charges ac- 
counting for the sign of the product gg-j- • Since the whole system must be "electrically neutral", at 
this point it is possible to develop the Debyc-Hiickel-type theory of hydrodynamic screening by analogy 
with that developed for Coulombic systems, e.g. see Ref.[41]. Nevertheless, below we choose another, 
more elegant pathway to arrive at the same conclusions. 



*Using dimensional analysis performed for Eq.(4.36) the result, Eq.(4.38), can be easily rewritten also in dimensionless 
form 
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Before doing so, we notice that there is an important difference between the double integral, Eq.(4.39), 
t 

and / [r(Tj)] drj present in the exponent in Eq.(4.31). While the double integral, Eq.(4.39), is 
° 

manifestly reparametrization invariant, the diffusion integral is not. This means that we can always 
reparametrize time in this diffusion integral so that the coefficient (4D)~^ can be made equal to any 
preassigned nonnegative integer. This was effectively done already when we introduced the dimensionless 
variables in Eq.(4.36). Such inequivalence between these two types of integrals can be eliminated if we 
replace this diffusion -type integral by that which is manifestly reparametrization- invariant. In such a 
case the total action is given by 

It should be noted that use of a symbol ^ instead of J in Eq.(4.40) is a delicate matter. In [33] we 
demonstrated that in the limit of long times (that is in the limit uj used in this work) all random 
walks are asymptotically closed (that is, the Brownian trajectory in this limit becomes very much the 
same as known for ring polymers)^. Since the result, Eq.(4.40), is manifestly reparametrization invariant, 
such a replacement is permissible. Additional explanations are given in Appendix B which we recommend 
to read only after reading of Section 5. 

The constant ttiq in Eq.(4.40) will be determined in the next section. The form of the action given 
by Eq.(4.40) is almost identical to that for the action for the superfluid liquid ^He as discussed in the 
book by Kleinert [42], page 300. From the same book, it also follows that the Ginzburg-Landau theory 
of superconductivity also can be recast in the same form. We said "almost identical to" meaning that in 
these two theories (of superfluidity and superconductivity) the self-interaction of vortices is also allowed 
so that if the above expression would represent the dual (vortex) description of colloidal suspension 
dynamics (e.g. see Appendix B), then the prime in the double summation above can be removed since 
the vortices are allowed to intersect with themselves. 

In the direct case, when the focus of attention is on particles, removal of the prime in the double 
summation in Eq.(4.40) would imply that the Oseen tensor is defined for particles hydrodynamically 
interacting with themselves. This assumption is not present in the original Doi-Edwards formulation, 
Ref. [10]. As we noticed already in Eq.(3.29), the diagonal part of the Oseen tensor is associated with self- 
diffusion. The question therefore arises: can this "almost equivalence" be converted into full equivalence? 
The main feature of superconductors is the existence of the Meissner (for hard spheres) and dual (for 
vortices) Meissner effect. In the present case such an effect is equivalent to the existence of hydrodynamic 
screening. Hence, to prove such an equivalence requires us to prove the existence of hydrodynamic 
screening for suspensions. Evidently, we cannot immediately use Eq.(4.40) for such a proof. Therefore, 
in the next subsection we use London-style arguments to arrive at the desired conclusion. 

4.4 London-style theory of hydrodynamic screening 

We begin our proof by taking into account the non-slip boundary condition, Eq.(2.27): 

v(r,t) = |=v(t). (2.27) 

Within the approximations made, we also have to impose the incompressibility requirement 

V-v(r,i)=0. (3.14) 

Because of this requirement, the current j = pv becomes j = riov with the density no being a constant. 
Since j is a vector, we can always represent it as 

j = aVi; (4.41) 
^Additional mathematical results on this property are discussed in Section 6.2. 
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with suitably chosen scalar ip and some proportionality constant a. To choose such a scalar we take into 
account that in the present case 

V-j=0 (4.42) 

implying 

VV = 0. (4.43) 
The vector j given by Eq.(4.41) is not uniquely defined. It will still obey the Eq.(4.42) if we write 

j = aV->p ± gA (4.44) 

for a vector A such that V • A =0. Evidently, a vector obeying Eq.(4.13) by construction possess this 
property. The choice of the sign "+" or "-" in the above equation can be determined based on the 
following arguments. Since j = nov and since no is constant, we can replace Eq.(4.44) by 

V = aVip ± gA (4.45) 

by suitably redefinig constants a and g. Next, we assume that v is a random variable so that on average 
< V >= thus implying 

< aVip >±g<A >=0. (4.46) 
This equation causes us to choose the sign After this, we can write for the correlator 

< V • V >= < V'^ • VV' > +g^ < A ■ A >= 2g^ < A ■ A > . (4.47) 

In view of our choice of A, the < A- A > correlator coincides with that given in the exponent of Eq. (4.24). 
Now wc take into account Eq. (4.20b) where, of course, we replace j by v so that using the dictionary, 
Eq.(4.16), we arrive at 

u) = ev (London equation) (4.48) 

supplemented with 

a; = V X V (Maxwell equation). (4-14) 

Such an identification becomes apparent because of the following arguments. Let us use Eq.(4.45) in 
Eq.(4.48) in order to obtain 

u} = e(aVV' - eA). (4.49) 

In this equation we replaced the constant g by e. Furthermore, since Eq.(4.49) formally looks like the 
Pick's first law, we can as well rewrite this result as 

u, = ei^Vi; - eA). (4.50) 

By applying to both sides of this equation the curl operator and taking into account Eq.(4.13), we obtain 

V X w = -e^v. (4.51) 
Taking into account the Maxwell's Eq.(4.14) and using it in Eq.(4.51) we obtain as well 

V^v = e^v. (4.52a) 

Equivalently, we obtain, 

V^A =e^A. (4.51b) 

Using Eq.s(4.47), (4.52a) and following the same steps as in the Appendix A of our previous work, Ref.[23], 
we obtain 

(v(r).v(0)) = ^exp(--^), (4.53) 
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( 



1 [D^ 



-1 



where = e 



-1 



i?o V 47rr7 



and the constant in Eq.(4.53) can be obtained from comparison 



between this equation and Eq.(2.37). The analogous result is also obtained for the < A • A > correlator. 

In accord with Eq.(2.44) we obtain the result of central importance — » oo when p ^ 0, imply- 
ing absence of screening in the infinite dilution limit. Our derivation explains the rationale behind the 
identification of Eq.s (4.14) and (4.48) with the Maxwell and London equations in the theory of super- 
conductivity, Ref.[25], pages 174, 175. Evidently, such an identification becomes possible only in view of 
the topological nature of the London equation, Eq.(4.48), coming from identification of Eq.(4.19) with 



5 Exotic superconductivity of colloidal suspensions 

5.1 General Remcirks 

In the previous section wc developed a theory of hydrodynamic screening following ideas of the London 
brothers, Ref.[24]. As is well known, their seminal work found its most notable application in the theory 
of ordinary superconductors [25]. At the same time, Eq.(4.40) was originally used in the theory of 
superfluid ^Hc. In the book by Klcincrt [42] it is shown that Eq.(4.40) can be rewritten in such a way 
that it will acquire the same form as used in the phenomenological Ginzburg-Landau (G-L) theory of 
superconductivity [25]. We would like to arrive at the same conclusions differently. In doing so we also 
would like to determine both the physical and mathematical meaning of the parameter mo which was left 
undetermined in Eq.(4.40). We shall develop our arguments mainly following the original G-L pathway. 

It should be said, though, that in the present case the connections with superconductivity are only in 
the structure of equations to be derived. The underlying physics is similar but not identical to that for 
superconductors. Indeed, in the case of superconductors one typically is talking about the supeconducting- 
to-normal transition controlled by temperature. Also, one is talking about the temperature-dependent 
"critical" magnetic field (the upper and the lower critical magnetic fields in the case of superconductors of 
the second kind) which destroys the superconductivity. In the present case of colloidal suspensions there 
is no explicit temperature dependence: the same phenomena can take place at any temperature at which 
the solvent is not frozen. If wc account for short range forces, then, of course, one can study a situation 
in which such a colloidal suspension is undergoing a temperature-controlled phase transition. Such a case 
requires a separate treatment and will not be considered in this work. In the present case the phase 
diagram can be qualitatively described as follows. The infinite dilution limit corresponds to the normal 
state. The regime of finite concentrations corresponds to a mixed state, typical for superconductors of the 
second kind, and the dramatic jump in viscosity discussed in the Introduction and in Section 2 corresponds 
to the transition to the "fully superconducting" state. Such a difference from the usual superconductors 
brings some new physics into play which may be useful, in other disciplines, e.g. in the high energy 
physics or turbulence, etc.^'^ 

5.2 G-L style derivation of equations of superconductivity for colloidal sus- 



We begin with the one of Maxwell's equations in its conventional form, e.g. as given in Ref.[25], page 181, 



(4.20a). 



pensions 




(5.1) 



In the G-L theory we have for the current j the following result: 



j = -7i-(<^ - ^w^*) — i^r A. 



(5.2) 



lOE 

.g. see Section 6. 
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Both equations can be obtained by minimization of the following (truncated) G-L functional^ ^ 



Stt 4m 



(5.3) 



with respect to A. Substitution of the ansatz ^p 



■ exp{itl)) into Eq.(5.2) leads to the current 



eh , 2e ^ , 



(5.4a) 



to be compared with our Eq.(4.50). Evidently, this result is equivalent to the postulated London equation 

for superconductors 



V X j = 



en. 



B. 



(5.4b) 



At the same time, a comparison of Eq. (5.4a) with Eq. (4.50) leads to the following chain of identifications: 

elx e? fi^ eng e^Tt 

ng ^ eDo and — ng ^ e^. Consequently, we obtain as well: - — ^ I?n, — > e; > e^,c 



2m 



2 m, 



h 



mc 

1 2e e 2e 

47r ^ — , — <=^ -— — . Using these identifications, we can rewrite the functional J^[A,(f] as follows 
27r n Do ris 



J-[A,^] = ^ J d^iVxAf+Do 



(V-i^A)^ 



(5.5) 



In the traditional setting, the superconducting density is determined from the full G-L functional 



^[A,^] = / £r{ 



(V X A)" 



Stt 



+ 



4m 



(5.6) 



e.g. by minimization with respect to (p*. In fact, to obtain rig it is formally sufficient to treat only the 
case when A = 0. Indeed, under this condition we obtain 



(5.7) 



which has a nontrivial solution only for a < 0. In this case we get Ug = provided that 6 > 0, as usual. 



■ exp(i'^) in Eq.(5.6) then, the polynomial 



If we use this result back in Eq.(5.6), that is we use (p^ = 

(in ip) part of the functional becomes a constant. This constant is divergent when the volume of the 
system goes to infinity. To prevent this from happening another constant term is typically added to the 
functional ^[A, tp] so that it acquires the following canonical form 



,3„,(VxA)^ ^ 

Stt 4m 



2ie ^ , 



(5.8) 



Then, when (p^. = 
in this limit we require 



cxp(i'0), the polynomial (in i^) part of the functional vanishes and, accordingly. 



/ 



2ie ^ , 



as well. This leads us to the equation 



or to 



— — Vip^ = A 
i2e 



he 
2i 



V^ = A. 



(5.9) 

(5.10a) 
(5.10b) 



^^This truncation is known in literature as the "London limit" . 
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This equation coincides (on average) with the previously obtained Eq.(4.46) (with redefinitions described 
above) as required and will be treated further in Section 5.4. 

It should be noted though that originally, in London's theory, Ref.[24], the rig was left as an adjustable 
parameter and, hence, microscopically undefined. This is important in our case since the phenomenon 
of supercoductivity can be looked upon (as in thermodynamics) without any reference to spontaneous 
symmetry breaking, Higgs effect, etc. At the level of G-L theory, the London equations are reproduced 
with help of the truncated G-L functional. Hence, in principle, in the present case use of the truncated 
fimctional, Eq.(5.5), is also sufficient. At the macroscopic mean field level the presence of polynomial 
terms in the fuU-G-L functional, Eq.s(5.6) and (5.8) seems somewhat artificial. They do not reveal their 
microscopic origin and are introduced just to fit the data. We would like to use some known facts from 
the path integral treatments of superconductivity/ superfluidity in order to reveal their physical meaning. 
Such information is also useful for development of the hydrodynamic theory of colloidal suspensions. 



5.3 Path integrals associated with the G-L functional 



In view of Eq.(4.40), we begin our discussion with the simplest case of the path integral for a single 
" relativistic" scalar particle. 

Following Polyakov, Ref.[43], the Euclidean version of propagator for such a (Klein-Gordon) particle 
is given by 

1 

G{x,x') = J (1^) exp(-mo J rfr^i^), (5.11a) 



where in the most general case 

x2(r)=5,.(x)^^. (5.11b) 

This propagator is of interest in string theory since it represents a reduced form of the propagator for the 
bosonic string. As in the case of a string, the action of this path integral is manifesttly reparametrization- 
invariant, i.e. invariant under changes of the type x(t) — > x(/(r)) (with /(r) being some nonnegative 
monotonically increasing function) . The path integral measure is designed to absorb this redundancy. The 
full account of this absorption is cumbersome. Because of this, instead of copying Polyakov's treatment of 
such a path integral, we shall adopt a simplified treatment allowing us to recover Polyakov's final results. 
We begin with an obvious well-known identity 

x(0)=0 

This identity is used below as follows. Consider the propagator for the Klein-Gordon (K-G) field given 
By employing the identity 

oo 

— = J dxexp{—ax) (5-14) 
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Eq.(5.13) can be rewritten as follows 



/f d k 
dt exp(-fm^) j exp(ik • x-tk^) 





oo 



= 4 / dtexx>{—£tm^) f ^ , exp(ik • x— gfk^) 

°° - 2 



= ^JdteM-tm') J S)[x(r)]exp{-^ydr(^j } (5.15) 

x(0)=0 

where we used the identity, Eq.(5.12). to obtain the last line and introduced an arbitrary nonnegative 
parameter £ for comparison with results by Polyakov. Specifically, using page 163 of the book by Polyakov 
(and comparing our 3rd line above with the 3rd line of his Eq.(9.63)) we can make the following identifi- 
cations: £ ^ e, w? ^ ji. Since, according to Polyakov, n = e~^(mo ^) with c being some constant, 

we obtain: mo = fm^ H — r=. That is, the physical mass entering the K-G equation is obtained as the 
limit of the expression {e — > 0) 

m= lim e~2 (mo — mcr)^- (5.16) 

mo — ^rricr 

Clearly, such an expression is nonnegative by construction. Prom the last line of Eq.(5.15) it follows that 

the propagator for the K-G field is just the direct Laplace transform of the nonrelativistic "diffusion" 
propagator, Eq.(5.12), with the Laplace variable m playing a role of a mass for such a field. In the 
Euclidean version of the K-G propagator this mass cannot be negative since in such a case the identity 
Eq.(5.14) cannot be used so that the connection between the nonrelativistic and the K-G propagators 
is lost. However, Eq.(5.2) seemingly is for the quantum current while the propagator in Eq.(5.12) is 
describing Brownian motion, not quantum diffusion. To fix the problem we have to replace time t in 
Eq.(5.12) by it and, accordingly, to make changes in Eq.(5.15). This then converts the Laplace transform 
into the Fourier transform, provided that the nonrelativistic propagator describes the retarded Green's 
function. To use the full strength of the apparatus of quantum field theory one needs to use the causal 
Green's functions. This is required by the rclativistic covariancc treating space and time coordinates 
on the same footing. Once all of these requirements are met, it becomes possible to treat the case of a 
negative mass. 

It should be emphasized at this point that the London-style derivation given in the previous section 
formally does not require such quantum mechanical analogy. Because of this, the following problem 
emerges: is it possible to reproduce the functional integral defined by Eq.(4.31) using the truncated 
G-L functional for superconductivity in the exponent of the associated path integral? We would like to 
provide an affirmative answer to this question now. 

We begin with the partition function Z for the two-component scalar K-G-type field 

lnZ= -In[det(-V2-Fm2)] (5.17) 

Since 

ln[det(-V^ + m^)] = tr [ln(-V^ + m^)] (5.18) 

and 

tr [ln(-V' + TO^)] = / ln(k2 + m^), (5.19) 
J (27r) 
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we can use the results of our previous work, Ref.[44], for evaluation of the last expression. Thus, we 
obtain. 



In 



(-V" + m2) 



(-v^) 



^ ' 

2 2 

m m 
^ ' 

oo x(t)=0 t 2 

= j dt j dyeM-ty) J 2D[x(r)]exp{-i J dr (^^^ } 

x(0)=0 

oo x(t)=0 t 2 

= J^il-eM-mh)) J S)[x(r)]exp{-iy'dT(^^) }. (5.20) 

x(0)=0 

x(t)=0 

Following the usual practice, wc shall write f instead of / in the path integral and consider a formal 

x(0)=0 

(that is diverging!) expression for the free energy J^q 



dt 



exp(-J-o) = In^o = -ln[det(-V^ + to2)] = / — exp(-m^f) f S)[x(r)] exp{-- j dr ( — ] } (5.21) 



1 



by keeping in mind that this result makes sense mathematically only when the same expression with 
TO^ = is subtracted from it as required by Eq.(5.20). Inclusion of the electromagnetic field into this 
scheme can be readily accomplished now. For this purpose we replace the V operator by its covariant 
derivative:V — > D = V— ieA (we put Dq = 1 in view of developments presented in Eq.(5.15)) Using D 
instead of V in Eq.(5.20) we have to evaluate the following path integral^^ 



[dct{-D^ + m^)]-^ = J D[ip,ip] cxp(-i y d^r{ip{-D'^ +m'^)ip}). 



(5.22) 



For A = we did this already while for A 7^ we can treat terms containing A as perturbation. We can 
do the same for the path integral in Eq.(4.32). This is easy to understand if we realize that 

00 00 * 2 

J dtexp{-mH)l[A;t] |a=o= J y exp(-m2t) ^ S)[x(r)] exp{- J j dr (^^^ } (5.23) 



1 fdxV 



Therefore, the final answer reads as follows 

oc t 

exp {-J^) = InZ = - ln[det(-D2 + m^)] = J ^ exp(-m^t) ^ S)[x(r)] exp{- j dT[ 


00 * 2 

= J y exp(-m2i)^2)[x(T)]exp{-^ J dr (^^^ } exp{-i^ ^dr • A[x(r)]}. 



4\dr) 



(5.24) 





00 



This result demonstrates that applying the operator / ^ exp(— m^f) to I[A; t] defined in Eq.(4.32) makes 



it equivalent to the "matter" part of truncated G-L functional for superconductivity as needed. This raises 



^^For m? = this is just part of the truncated G-L functional. 



28 



a question about comparison of the full G-L functional with the "diffusion" path integrals of Section 4. 
Evidently, this can be done only if in the original diffusion Eq.(3.30) we do not discard the potential U. If 
we do not discard the potential and if, instead, we ignore the hydrodynamic interactions completely, we 
would end up with the following path integral for interacting Brownian particles in the canonical ensemble 



E= l[V[^{ri)] exp{-— ^ / dr, f ^ j ^T. ^tJ dr,-F[x(r,) - x(r,-)]}. (5.25) 

It is essential that this expression does not contain self-interactions typical for problems involving polymer 
chains with excluded volume-type interactions. The situation here resembles that encountered when, 
following Doi and Edwards, Ref.[10], we redefined the Oseen tensor in Eq.s (3.29 ) and (3.30) so that it 
acquired the diagonal part. In the present case we must require the diagonal part to be zero at the end 
of calculations. These results, correct for colloidal particles, may become incorrect in the present case 
for the following reason. Prom looking either at Eq.(5.24) or Eq.(4.39), we recognize that in these cases 
we are dealing with assemblies of loops (vortices) which are in one-to one correspondence with diffusing 
particles. While this topic is studied in detail in the next subsection and Appendix B, here we notice 
that if Eq.(5.25) is written for such loops, then the excluded volume requirement becomes essential, even 
for a single loop. Indeed, the existence of such a loop is possible only if the field A associated with those 
loops is uniquely defined. This is possible only if the loop contour does not have self-interactions. This 
is the origin of the excluded volume constraint requirement. With this restriction imposed, we introduce 
a density 

N * 
'=1 

so that the binary potential in Eq.(5.25) can be written as 

N ^ ^ 

^ E / / dTiV[^{Ti) - x(r,)] = \jdvj dv'p{v)V[v - r>(r'). (5.27) 

^'i 

Then, using the Hubbard-Stratonovich (H-S) identity we obtain, 

exp(-iy" dv j dv' p{y)V[v-v'\p{y')) j D[iP{r)]exp{-^ j dr j (ir'V'(r)y-i[r-r']V'(r')) exp(i J drtl;{r)p{r)) 

(5.28) 

with or being a normalization constant (bringing the above identity to the statement 1 = 1 for p = 0). 
Use of this result in Eq.(5.25) in which self exclusion is allowed converts this partition function into the 
following form (written for the loop ensemble) 

N 

S = <ny £>[V(r)]exp(-i J dr J rfrV(r)y-i[r - r']V(r')) Gi(0; i | V) (5.29) 

where 

Gi{0;t I V) = y S)[x(r,)]exp{- J dnl^ (^^^ + ie^[x(TO]]}. (5.30) 



In the case of polymers, typically, one uses the delta function-type potential for description of the inter- 
actions. This observation is helpful in the present case as well because of the following. Consider the G-L 
functional, Eq.(5.6), and use the H-S identity for the interaction term 

exp(-^[|^|2]2) = nj i?[^(r)]exp(-l J dr J (irV(r)^(r')) exp(z J dr^{r)) \^\^). (5.31) 
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This allows us to replace the determinant, Eq.(5.22), by the following (more general) determinant 

[det(-D2 + TO^ + i^p)]-^ = j D[ip, ip] exp(-i J d^r{ip{-D^ +m^ + #[r])(p}) (5.32) 
which, in view of Eq. (5.24), can be equivalently rewritten as 

exp (-J^) = InZ = - ln[det(-D2 + + #)] = mj D[ip{r)] exp(-^ / / dr'i'{r)ip{r')) 

oo * 2 

X J y exp(-m^t) ^2)[x(t)] exp{-i J } exp{-ie j) dx ■ A[x{t)]} + i j) dT^p[x{T)]} 



oo t ^ 

= j j^M-mH) j^[x{T)\ex^{-^j dT(^ 





X exr 



} 



sp{-iy j dx ■ A[x(t)] -\jdTj> dT'5{x{T) - x(t')}. (5.33) 
Alternatively, this result can be rewritten as a grand canonical ensemble of selfavoiding loops 



dx 
dri 



} 



Z[A;b]-l = ^-^n[/|^exp(-m2f0/D[x(r0]exp{-iX:/rf^< 
Ti=l ■ z=i ;^ ' 1=1 

X exp{-i-^ / dx • A[x(r;)]} " | i dn i dT'^6{x{Ti) - x{t'^)}. (5.34) 

■' •' l,m=l-' ■' 

This result is useful to compare with Eq.(4.31). Prom such a comparison it is evident that Eq.(5.34) 

is compatible with that obtained previously. It accounts for the effects of non hydrodynamic-typc in- 
teractions which can be incorporated, in principle, in the diffusion Eq.(3.30) in which the potential U 
must be specified. Clearly, the use of path integrals makes such a task much simpler. However, even 
though the above derivation is intuitively appealing, strictly speaking, it cannot be used for a number of 
reasons. Unlike the G-L functional, Eq.(5.8), which is convenient for studying of topological and nonper- 
turbative effects, Eq.(5.34) makes sense only in perturbative calculations. This means that phenomena 
such as screening (caused by the Higgs effect) cannot be captured with such a formalism alone. These 
observations explain why screening effects were found in solutions of polymers but not in colloidal sus- 
pensions, Ref. [10,12]. Furthermore, Eq.(5.34) contains a mixture of reparametrization-invariant and non 
invariant terms. This is questionable mathematically. It would be more logical to have the entire action 
reparametrization-invariant. We study these issues in some detail in the next subsection. 

5.4 Repeirametrization-inveiriance and its consequences. London-style analy- 
sis 

Since path integrals mathematically can seldom be defined rigorously, we would like in this subsection to 
extend the analysis of Sections 4.2.-4.4 avoiding the use of path integrals. We start with a discussion of the 
result, Eq.(5.10b), which is the superconducting analog of Eq.(4.46) for suspensions. Since B = V x A, 
(or V = V X A in the case of suspensions) we conclude that Eq.(5.10b) causes the first term in the G-L 
functional Eq.(5.6) (or (5.5)) to vanish in the bulk. Nevertheless, it is perfectly permissible to write 



2e 

and to use Stokes' theorem 



^ (hdr-Vi}= (hdr-A. (5.35) 



dr-A = / / dS • (V X A) = / / dS • B =n^ (5.36) 



c 
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with n = 0,±1,±2,... An analogous result for suspensions reads as follows 

jdT-A= JJ dS- {V X A) = JJ d^-"^ (5.37a) 

c 

In view of Eq.s(4.19) and (4.20), this result leads also to 



e 



dS • w =n = — (/) dr • A (5.37b) 



c 



which is the same as Eq.(4.21). 

These results can be interpreted in a number of ways. For the sake of argument, we would like to 

explore the more established case of superconductivity first. Following Lund and Regge, Ref.[45], we 
suppose that the vector potential A can be presented as follows 

A(r) = A / 1^ (^\da=^i ^—-v{a)d<j (5.38) 



A-K J |r — r((7)| \da J An J |r — r(cr)| 

c c 

with the appropriately chosen constant k. This result easily follows from Eq.(4.15) imdcr the assumption 
that (Z)(r) — k f dav{(T)6{r — r{a)) (which is the same as our Eq.(4.20b)). Substitution of this result into 

c 

Eq.(5.36) produces 

. A = A / / dada',^^14-^ = n'^. (5.39) 
47r// |r(o-)-r(c7')| 2e ^ ' 

The obtained result allows us to determine the constant k. To do so we need to demonstrate that the 
above double integral is a linking number, e.g. see Eq.(4.6). The proof of this result depends upon 

correctness of the following statement 

dr ■ A = dr • B dr • V X A = 

c c 



47r/ ^ ' J ^ ' |r(o-)-r(c7')|'' 



Ci C2 



^ttJJ |r(<7) - r(a')| 

= Kklk{l,2) (5.40) 

with linking number lk{l, 2) defined in Eq.(4.6). If the above result is correct and the constant >f can be 
found then, the constant k can be determined from Eq.(5.39). Hence, the task lies in demonstrating that 

the nonzero constant x does exist. To do so we shall use the standard London analysis. Thus, we write 

47r 

V X B = — j (Maxwell equation) (5.1) 
c 

and 

Vxj = ^B (London equation). (5.4b) 

Since, B = V x A, and V • B = V • A =0, we have xB = A so that we obtain 

V X A =x— j (5.41) 
c 



and, from here 



-V^A = V X V X A =x— (V X j) = -X r^B = -x^ r^A (5.42) 
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which is the familiar screening-type equation, e.g. sec Eq. (4.51b). Since, in the conventional setting the 
penetration depth 6 is known to be (5 = ( ^ ) , we can chose implying that k — —. The 



\ mc^ / 2e 
choice >f=l does not mean of course that the constant k is dimensionless. Because of this, we obtain 

(5-43) 

in accord with Eq.(4.21). Next, if wc take into account screening effects, the conclusions we've reached 
will remain the same due to reparametrization invariance of both sides of Eq.(5.43). Indeed, consider one 
loop, say Ci, going from -oo to +00 in the z-direction. If we compactify by adding one point at infinity 
so that T{? becomes S"^, then such a loop will be closed. Another loop can stay mainly in the x-y plane so 
that the linking number becomes the winding number, e.g. see Ref.[46], page 134. Under these conditions 
the screening factor exp{—-) ^^under the integral of the left hand side of Eq.(5.43) is unimportant since 

we can always arrange our windings in such a way that r <i; 5 for any preassigned nonzero 8 so that the 
screening factor becomes unimportant. 

The above analysis can be extended to the case of colloidal suspensions in view of the results of Sections 
4.2 and 4.4. implying that in both superconductivity and colloidal suspensions the phase transition is 
topological in nature (e.g. in the colloidal case Eq.(4.39) is a topological invariant to be considered in 
the next subsection). Evidently, such a conclusion cannot be reached by perturbatively calculating the 
Green's function in Eq.(4.37). 

In Section 5.1 we discussed similarities and differences between superconductors and colloidal sus- 
pensions. It is appropriate now to add a few additional details to the emerging picture. In the case of 
superconductivity correctness of the topological picture depends upon the existence of nontrivial solu- 
tions of Eq.(5.42). These arc possible only when the parameter n,, is nonzero. When it becomes zero the 
above picture breaks down. In the case of suspensions the role of the parameter >c~^ is played by the 
density-dependent parameter e. This can be easily seen if we take into account that dimensional analysis 
requires us to replace Eq.(5.38) by 



c 



so that by employing Eq. (5.37b) we obtain, 



C\ C2 

as expected. 



^ i i dcrdcr' , ^j^'^} ' ^^1, =n = Ml, 2) (5.45) 
An J J |r(o-)-r(c7')| 



5.5 Bose-Einstein-type transition in a system of linked loops 

In the Introduction we noted that Chorin, Ref.[22], conjectured that the superfluid-to normal transition 

in '*IIc is associated with vortices causing a sharp increase in viscosity. In this subsection we wold like 
to demonstrate that, at least for colloidal suspensions, his conjecture is correct: the sharp increase in 
viscosity is associated with the lambda-type transition. Instead of treating this problem in full generality, 
i.e. for the nonideal Bose gas, we simplify matters and consider a Bosc condensation type transition 
typical for the ideal Bose gas. It should be noted though that our simplified treatment is motivated only 
by the fact that it happens to be sufficient for comparison with experimental data. In other cases, such 
a restriction can be lifted. 

^^Emergence of such a screening factor can be easily understood if we replace Eq.(4.15) by Eq.(5.42) with the right hand 
side given by w(r) = k § v{a)5{T — r(cr)) in accord with Eq.(5.38). 

C 
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To develop such a theory we use the information obtained in the previous subsection augmented by 
some additional facts needed for completion of our task. In particular, we are interested in the expression 
for the kinetic energy. Up to a constant it is given by 

E=^ j Sr{V X Af (5.46) 

and is manifestly nonnegative. Using known facts from vector analysis this expression can be rewritten 
as follows 

E=^ j d^r{V X A) • V =^ j rf^r[A • ib+div[A x v]] =^ j d^rA ■ w (5.47) 
In view of Eq.(5.38) we can rewrite this result as 

E^^m dada' (5.48) 



2 J J \v{a) - v{a')\ 

Ci C2 

to be compared with Eq.(5.45). Using such a comparison we arrive at an apparent contradiction: while 
an expression for E should be nonnegative, the linking number lk{l, 2) can be both positive or negative. 
If we make the replacement r — r in Eq.(5.48) nothing changes but if we do the same for lk{l,2) it 
changes the sign. Thus, if we want to use ZA:(1,2) in Eq.(5.48) we have to use |^fc(l,2)|. This number 
was introduced by Arnold and is known in literature as entanglement complexity^'^ . Evidently, in view 
of this remark, n in Eq.(5.45) can be only nonnegative. If wc require our system to be invariant with 
respect to rotations of the coordinate frame, Eq.(4.39) should be rewritten according to the procedure 
developed in our work, Ref.[49]. This means, that we introduce a set of linking numbers: ni, n2, rzi, ... 
so that for a given n^^ , the set of ^n{n — 1) = N possible linking numbers can be characterized by the 
total linking number L, i.e. we have 

^ni = L. (5.49) 

i=l 

This result can be rewritten alternatively as follows. Let Ci be the number of links with linking number 
1, C2 the number of links with linking number 2 and so on. Then, we obtain 

L 

^iCi=L. (5.50) 
1=1 

Furthermore, we also must require 

L 

Y^Ci = N (5.51) 

i=l 

Define the Stirling-type number S{L, N) via the following generating function^^ 

L 

S{L, N)x^ = x{x + l)---{x + L-l). (5.52) 

N=0 

Set in this definition a; = 1. This then allows us to introduce the probability p(L, A^) = S{L,N)/Ll The 
number S{L, N) can be easily obtained^'' with the result given by 

^ 7-1 

^(^'^) = n^- (5.53) 

i—l * 



^^For more dcatails about this number and its many applications can be found in our works, Refs.[47,4S 
l^E.g. see Eq.(4.4). 

i^The true Stirling number of the first kind S(L, N) is defined as follows; S{L, N) := (-1)^-^5(L, N). 
i^'E.g. see Ref. [49]. 
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With thus obtained results, we are now ready to return to Eq.(4.39) in which we make a rescahng: 
r(T) — ^ i?of(r), with f(r) being dimensionless. After which, using Eq.(5.45) we can rewrite Eq.(4.39) as 
follows 

< WiL) >T= exp(-^L) (5.54) 

Evidently, the numerical factor of 3 in the exponent is non-essential and can be safely dropped upon 
rescaling of L. To use this expression we combine it with Eq.(5.34) in which we have to make some 
adjustments following Feynman, Ref.[50], pages 62-64. On these pages Feynman discusses a partition 
function for the ideal Bose gas written in the path integral form. We would like to rewrite his result in 
the notation of our paper. For this purpose we use Eq.(4.30) in which the path integral is written for a 
loop and is in discrete form. We obtain, 

^^""^ = (i^) ^ /n^'r,exp{^^[(ri-r2)^ + (r2-r3)^ + ... + (r._i-r.)^ + (r.-rin 

with V ~ J d^ri. Under such circumstances the Brownian ring is made out of v links (segments) so that we 
can identify its length with v. In the present case each such ring is linked with another ring thus forming 
a link with a linking number iCt, i = 0, 1,2, ... Since the linking number is independent of the lengths 
of rings from which it is made, we can take advantage of this fact by identifying the index i with v. By 
combining Eq.s (5.50)-(5.55) and repeating the same steps as given in Feynman's lectures we assemble 
the following dimensionless grand partition function 

Ci,...,C„.... f " ' Ci,...,C,,.... f 

Here the "chemical" potential z = exp{—^). Taking the logarithm of both sides of the above equation 
we obtain the partition function for the ideal Bose gas. Written per unit volume it reads 

3 

^=-(^)^C.M.). (5.57) 
In this expression (^^ (z) is Riemann's zeta function 



C«W = E;;^- (5-58) 



This function is well defined for 2; < 1, i.e. for — < 00 and is divergent for 2; > 1, thus indicating a Bose 

Vo 

condensation whose onset is determined by the value z = 1 (i.e. rj = 00) for which C5/2(l) = 1-341. If we 
follow standard treatments, then we obtain for the critical density 



4ttDo 



2.612. (5.59) 



In view of Eq.(5.55), the obtained result for density has the correct dimensionality. Prom here the critical 
volume fraction is: (fi^ = p^^ttRq. The number 2.612 is just the value of C3/2(l)- This means, that we can 
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write in the general case 

3 



thus giving us the equation 



^ = l_k44. (5.61) 

Pc C3/2(l) ^ ' 

In the book by London, Ref.[51], we found the following expansion for C^i2{z) in the vicinity of 2; = 1 
< 1) : 

(3/2(2) = -3.545a5 + 2.612 + 1.460a - 0.1040^ + (5.62) 
where a = — Inz. Use of this result in Eq.(5.61) produces the following result: 

^(^)'"-^^ 

in accord with scaling predictions by Brady, Ref. [19], and Bicerano et aLRef.[20]. It should be noted 
though that in view of Eq.(5.54) the actual value of the constant prefactor in Eq.(5.63) is quite arbitrary 
and can be adjusted with help of experimental data. For instance, by making this prefactor of order 
unity, Bicerano et all obtained a very good agreement with experimental data in the whole range of 
concentrations, e.g. see Ref. [20], Fig. 4. 



6 Discussion and outlook 
6.1 General comments 

With the exception of the work by De Genncs [52] on phase transition in smcctics A, the superconductiv- 
ity and superfluidity phenomena are typically associated with the domain of low temperature condensed 
matter physics^^.This fact remains true even with account of cuprate superconductors, Ref. [54]. The 

results obtained in this work cause us to look at these phenomena differently. For instance, the previously 
mentioned relation Cj{v) = k § w{a)5{Y — v{ij)) used in the work by Lund and Rcggc, Ref. [45], for fluids, 

c 

coincides with our Eq. (4.20b) for colloids. The work of Lund and Regge is based on previous work by 
Rasetti and Regge, Ref. [55], on superfluid He and, therefore, their results are apparently valid only in 
the domain of low temperatures. This conclusion is incorrect however as shown in the scries of papers 
by Berdichevsky, Ref.s [56,57]. Any ideal (that is Euler-type) incompressible fluid can be treated this 
way. Furthermore, as results by Chorin, Ref. [22], indicate, the same methods should be applicable for 
description of the onset of fluid/gas turbulence. In our work the fluid is manifestly nonideal. Nevertheless, 
in the long time (zero frequencies) limit it can still be treated as if it is ideal. 

The most spectacular departure from traditional view on the results by Lund and Regge was recently 
made in a series of papers by Schicf and collaborators, Rcfs. [58,59]. The latest results elaborating on 
his work can be found in Ref. [60] . Schief demonstrated that the results of Lund and Regge work well in 
the case of magnetohydrodynamics, that is, ultimately in the plasma installations designed for controlled 
thermonuclear synthesis. 

The basic underlying physics of all these phenomena can be summarized as follows. In every system 
which supports knotted structures, the existence of a decoupling of topological properties from the confor- 
mational (statistical) properties of flux tubes from which these knots/links are made should be possible. 
Since this statement is not restricted to a simple Abelian C-S field theory describing knots/links existing 
in G-L theory, in full generality the theory should include the G-L theory as a special case (as demon- 
strated above). Accordingly, the minimization of the corresponding truncated G-L functional may or 
may not lead to London-type equations. We would like to illustrate these general statements by specific 
examples. This is accomplished below. 

^^Lately, however, these ideas have began to be popular in color supercoductivity dealing with quark matter [53]. 
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6.2 Helicity and force-free fields imply knoting and linking but not nesess- 
SEirily superconductivity via London mechanism 

The concept of helicity has its origin in theory of neutrino, Ref.[61]. An expression a ■ p/ |p| is called 

helicity. Here cr • p = cr^p^ + fyPy + cr^Pz , and pi and di , i = 1 — 3, are being respectively the components 
of the momentum and Pauli matrices. The eigenvalue equation 

[o- • p/ |p|] * = A* (6.1) 

produces eigenvalues A which can be only ±1. Moffat, Rcf.[62], designed a classical analog of the helicity 
operator. He proposed to use the product v • V x v = v • d; for this classical analog. In it, as before, e.g. 
see Eq.(4.14), the vorticity field w is used. Moffat constructed an integral (over the volume M) 

V • u>dV (6.2) 

M 

along with two other integrals: the kinematic kinetic energy 

2T 



-E 

P 

M 



j v'^dV (6.3) 



and the rotational kinetic energy 



Then, he used the Schwarz inequality 



j u^dV. (6.4) 



O = y Gj^dV. 

M 



P 

I'^ <En or Q > — (6.5) 

E 

in order to demonstrate that the equality is achieved only if w = q;v where a is a constant. Since this 
requirement coincides exactly with our Eq. (4.20b), it is of interest to sudy this condition further. In 
particular, under this condition we obtain al = E which would coincide with our Eq.(5.43) (sec also 
5.48)) should / be associated with the linking number. Fortunately, this is indeed the case. The proof 
was given by Arnold and is outlined in Ref.[63], pages 141-146. In view of its physical significance, we 
would like to discuss it in some detail. 

Before doing so, we notice that the condition us = cxv is known in literature as the force-free condition 
for the following reason. In electrodynamics, the motion of an electron in a magnetic field is given by (in 
the system of units in which m = c = e = 1) 

^ = V X B (6.6a) 

dt ^ ' 

while the use of the Maxwell's equation, our Eq.(4.10), produces as well 

V = V X B = aB (6.6b) 

Using previously established equivalence v ^ B and substitution of Eq.(6.6b) into Eq.(6.6a) explains 
why the force-free condition is given by (Z) = ctv. This equation can be looked upon as an eigenvalue 
equation for the operator V x (• • •). From this point of view the force-free equation is totally analogous 
to its quantum counterpart, Eq.(6.1). Details can be found in Rcf.[64]. 

Going back to Arnold's proof, we note that according to Moffatt, Ref.[62], page 119, 

j = y,, .-.^T._ 1 /■ /■ Ri2 • P(l) X (2^(2)] 
V y(i)y(2) 
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Clearly, if as is done by Moffatt and others in physics literature (e.g. Lund and Regge, etc.), we assume 
that the vector potential A can be given in the form of Eq.(5.44), then / indeed becomes the linking 
number, Eq.(4.6). If, however, we do not make such an assumption, then much more sophisticated 
methods are required for the proof of this result. Use of these methods is not of academic interest only, 
as we would like to explain now. According to Kozlov. Ref.[65],thc force- free case u) = av belongs to the 
category of so called vortex motion in the weak sense. There are many other vortex motions for which 
V X V X V 7^ 0. These are vortex motions in the strong sense. Evidently, any relation with superconductivity 
or superfluidity (which is actually only hinted at this stage in view of results obtained in previous sections) 
is lost in this (strong) case. But even with the vorticity present in the weak sense this connection is not 
immediately clear. This is so because of multitude of solutions of the force-free equation as discussed, for 
example, in Rcfs.[66, 67]. Wc would like to discuss only those solutions which arc suitable for use in 
Arnold's theorem. These solutions can be obtained as follows. Taking the curl of the equation 

V X B = aB, (6.8) 

provided that V • B = 0, produces 

(V^ + a^)B = 0, (6.9) 

to be compared with our result, Eq. (4.52a). Unlike our case, which is motivated by analogies with super- 
conductivity and superfluidity, in the present case there are many solutions of this equation. We choose 
only the solution which illustrates the theorem by Arnold. It is given by v = {Asinz + Ccosy, Bsinx + 
Acosz, Csiny + Bcosx), where ABC ^ and A,B,C e R This solution is obtained for a = 1. 
Following Arnold, we introduce the asymptotic linking number A(a;i,a;2) via 

A(.„..)= lim -4- ? JdUdt. ^ • 3- . (6.10a) 

^ TuT,^^AnT,T2j J ||xi(ti)-X2(t2)f ^ ' 

II V / \ /II 

The theorem proven by Arnold states that if the motion described by trajectories xi(ti) and X2(t2) is 
ergodic, then 

1 [ f Ri2-p(l)xw(2)] 1 



, , , ^ '-dV{l)dV{2) = / / A(xi,X2)dV{l)dV{2)^lk{l,2). (6.10b) 

47r 7 J R12 ^ J J 

V{1) V{2) V{1) V{2) 

That is the function A(xi,a;2) on ergodic trajectories is almost everywhere constant. This theorem as 
such does not imply that this constant is an integer. For us it is important to realize that both Eq. (4.52a) 
and Eq.(6.9) can produce trajectories minimizing the Schwarz inequality thus leading to the condition 
al = E with I being either linking (in the case of suspensions) or self-linking number (depending upon 
the problem in question) or a conbination of both. Because both Eq. (4.52a) and (6.9) cause formation 
of links, the choice between them should be made on a case-by-case basis. In particular, existence of the 
Messner effect in superconductors leaves us with no freedom of choice between these two equations. In 
the case of magnctohydrodynamics/plasma physics the situation is less obvious. In the next subsection 
we shall argue in favour of superconducting/ superfluid choice between these equations. To our knowledge, 
such a choice was left unused in plasma physics literature. 

6.3 Ideal magnetohydrodynamics and superfluidity /superconductivity 

In order to discuss the work by Schief, Ref.[58], we would like to remind to our readers of some facts 
from the work by Lund and Regge (originally meant to describe superfluid '^He) since these fact nicely 
supplement those presented in previous sections. We already mentioned that Berdichevsky adopted these 
results for normal fluids, including those which are turbulent. Lund and Regge assumed that the vortex 
has a finite thickness so that the non-slip boundary condition, Eq.(2.27), should be now amended to 
account for finite thickness. The amended equation is given by 

dxi dxi da 
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where a parametrizes the coordinate along the vortex line. Eq.(5.38) taken from work by Lund and Regge 
then implies: 

%-fc^(^-^.) = 0. (6.12) 

This equation is treated as an equation of motion by Lund and Regge obtained with help of the following 
Lagrangian 

dxj dxk j_ p 



da dt 2 
c V 

Since the zero thickness limit of the action for this Lagrangian is given by our Eq.(4.22), which upon 
integration of the A-field leads to the result, Eq.(4.24), the same can be done in the present case and, 
accordingly, by analogy with the action, Eq.(4.22), which was extended, e.g. see Eq.(4.40), in the present 
case it can be extended as well so that the final result for the action of the Nambu-Goto bosonic string 
interacting with electromagnetic-type field reads (using the same signature of space-time as used in 
Ref.[45]) 

dx^ dx^ , , 1 
'da~dt' 



da-!- v^d^y. (6.13) 



S = -m I dadT^ + f I A^.^^dadT - - / F^dvol (6.14) 



with 



dx" dx„\ f dx^" dxa\ f dx" dx„\ /^irN 



da da I \ Or dr J \ da dr 



and m and / being some coupling constants. The metric of the surface enclosing the vortex can be 
always brought to diagonal form by some conformal transformation^^. In such coordinates, variation of 
the action S produces the following set of equations 

m(^--^)x =fs>^-^^PF^^ (6 16a) 

^Ot^ da^' " ^ " dr da ^ ' 

and 

d^dri^^ ~ ^^)S^^){xia,T) - y) (6.16b) 

provided that d^A^" = 0. Since the last equation is just the wave equation with an external source, the 
equation of motion for the vortex is Eq.(6.16a). In such a form it was obtained in Ref.[58] describing 
vortices in ideal magnetohydrodynamics. Under some physically plausible condition it was reduced in the 
same reference to the equation of motion for the one-dimensional Heisenberg ferromagnet. This result 
will be discussed further below from a somewhat different perspective. 

It should be noted though that Eq.(6.16a) emerges in Ref.[58] under somewhat broader conditions 
than those allowed by the force-free equation. In view of the content of the next subsection, we would 
like to reproduce this, more general case, now. For this purpose, we recall that the Euler's equation for 
fluids can be written in the form, Ref.[29], 

d 

—d} = V X (v X w). (6.17) 

In the case when u) is time-independent, it is sufficient to require only that 

V X w = V$ (6.18) 

with $ being some (potential) scalar function. In the case of hydrodynamics the equation $ = const is 
the famous Bernoulli equation. Thus, the force-free condition in this case is equivalent to the Bernoulli 
condition/equation. In magnetohydrodynamics there is an analog of the Bernoulli equation as explained 
in Ref.[69]. So, again, the equation $ = const is equivalent to the force-free condition. In the case of 
magnetohydrodynamics the vortex Eq.(6.16.a) is obtained under the condition $ = const. Since Eq.(6.16a) 
describes the vortex filament, the helicity integral, Eq.(6.7), describes either linking, self-linking or both. 



i^For more details, please see Ref.[68]. 
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In the case of self-linking it is known, e.g. see Ref.s[48,63], that lk{l, 1) = Tw + Wr. Analytically, the 
writhe Wr term is expressible as in Eq.(4.6) but with Ci and C2 now representing the same closed curve. 
The need for Tw disappears if the closed curve can be considered to have zero thickness. More accurately, 
the closed curve should be a ribbon in order to have a nonzero Tw. This is explained in Ref.[63]. With 
the exception of Appendix C, in this work we have ignored such complications. 

6.4 Classical mechanics in the vortex formalism, inertial dynamics of nonrigid 
bodies and G-L theory of high temperature superconductors 



Euler's Eq.(6.17) can be rewritten in the equivalent form: 

d 



at 



V X w ■ 



(6.19) 



Following Kozlov, Rcf. [65], in the case of Hamiltonian mechanics it is convenient to consider a very similar 
(Lamb) equation given by 

d 

— u+ (V X u) • v = ~ V$, (6.20) 

in which the vector u is such that V • u = 0. It can be demonstrated that Hamiltonian dynamics is 
isomorphic to the dynamics described by the above Lamb equation, provided that we make the following 
identifications. Let S" be a manifold in phase space P = T*M admitting a single-valued projection onto 
a configurational space M. In canonical coordinates x and y this manifold is defined by the equation 



y = u(x, t). 



(6.21) 



It is not difficult to demonstrate that the manifold S" is an invariant manifold for a canonical Hamiltonian 
H{x,y,t) if and only if the field y satisfies the Lamb's Eq.(6.20) and that $(x, i) = if(x, y(x, t), t) is a 
function on M parametrized by time t in such a way that 



dH 



(6.22) 



and 



dH , 

y = -— ly= 



(6.23) 



9u 5u 

= u 1 • V. 

9x " dt dx 

Relevance of these results to our discussion can be seen when Eq.(6.23) is compared with Eq.(6.11) of 
Lund and Regge. This comparison shows their near equivak;nc;c!. In view of this, wc would like to exploit 
this equivalence further by employing it for analysis of the truncated G-L functional analogous to our 
Eq.(5.3) typically used for phenomenological description of high temperature superconductors [54]. In 
this case the functional ^[A, ip] should be replaced by 



^[A,<p] 



(V X A)' 



+ 



2ie ^ . 
(Vx - ^Ax)^ 



-I- 



4m II 



2ie ^ , 
(V||-^A||)^ 



} 



(6.24) 



with its components lying in the x-y (cuprate) plane and z-plane perpendicular to it. By varying this 
functional with respect to Aj_ and Ay separately we obtain respectively the following components for the 
Maxwell's equation 

47r 

V X Bi = — ji (i =_L and II ), (6.25) 



where 



J-L = - 



^ (^* Vi(^ - Vi(^* )- ](/.]' Ax. and j II = 



r * ^ ^ *\ I |2 A fa ORA 
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Prom here we obtain the phenomenological London-type equations 

Vxj_L= ^B_L and V X j||= ^-Bn. (6.27) 

m_LC " m||C " 

By combining Eq.s (6.25) and (6.27) and using results of our Sections 4.2. and 4.4 we can rewrite these 
equations in the following suggestive (London-type) form 

u;_L = e_LVj_ and W|| = e||V||. (6.28) 

This form allows us to make a connection with the inertial dynamics of a nonrigid body. Following Kozlov, 
Rcf. [65], wc consider the motion of a nonrigid body in which particles can move relative to each other due 
to internal forces. Let the inertia axes of the body be the axes of the moving frame. Let K be the angular 
momentum of the body relative to a fixed point and u> the angular velocity of the moving trihedron while 
the inertia matrix I is diag (/_l, I±, /||)^°. The angular momentum and the angular velocity are related by 

K = Iw -I- A, (6.29) 

where A =(A_l, A_l, Ay) is the gyroscopic torque originating from the motion of particles inside the body. 
Prom here we obtain the Euler equation 

K = w X K =0, (6.30) 

which is a simple consequence of Eq.(6.29). In view of Eq.s(4.45) and (4.48) wc can identify Eq.s (6.28) 
with (6.29) thus formally making Eq.s (6.29), of London type. The hydrodynamic analogy can be in fact 
extended so that the hydro dynamically looking Lamb-type equation can be easily obtained and analyzed. 
Details are given in Ref.[65], page 148. 



6.5 Dirac monopoles, dual Meissner effect, Abelian projection for QCD and 
string models 

At this point our readers may have already noticed the following. Lin our derivation of Eq.(5.63) we 
made screening effects seemingly disappear while the title of our work involves screening. 2. In Eq.(6.14) 
wc introduced the Nambu-Goto string normally used in hadron physics associated with non Abelian 
Yang-Mills (Y-M) gauge fields. Quantum chromodynamics (QCD) of hadrons and mesons is definitely 
not the same thing as scalar electrodynamics (that is G-L model) discussed in our work. 3. Variation 
of the action S in Eq.(6.14) leading to the string equation of motion, Eq.(6.16a), under some conditions 
reduces to the equation of motion for the Heinsenberg (anti) ferromagnetic chain, which indeed describes 
the motion of the vortex filaments [59]. From this reference it follows that such equation of motion, in 
principle, can be obtained quite independently from the Nambu-Goto string, QCD, etc. In this subsection 
we demonstrate that the above loose ends are in fact indicative of the very deep underlying mathematics 
and physics needed for a unified description of all of these phenomena. 

The formalism developed thus far in this work suffers from a kind of asymmetry. On one hand, we 
started with a solution of hard spheres and then we noticed that these hard speres in solution act as 
currents (if one is using the magnetic analogy). The famous Biot-Savart law in magnetostatics is causing 
two currents to be entangled with each other thus creating the Gauss linking number, Eq.(4.8). Thus, it 
appears that in solution two particles (currents) are always linked (entangled) with each other. That this 
is indeed the case was noticed long ago as mentioned in the Introduction, e.g. sec Ref.[ll]. We can treat 
the vortices causing such linkages as independent objects. This is reflected in the fact that we introduced 
the vorticity (D (r) a.s u — k f da-v{(T)S{r ~ r{a)), e.g read comments after Eq. (5.38). In view of our major 

c 

equation a> (r) = ev, we can think either about the velocity (or vorticity) of a particular hard sphere or 
about the velocity of a particular vortex. Because of this, it is possible to treat both particles and vortices 

^"For the Scike of compaj-ison with superconductors, we deliberately choose the matrix in such form. 
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on the same footing. In such a picture (sketched in Appendix B) one can either eUminate vortices and 
think about effective interactions between hard spheres or vice versa. In this sense we can talk about the 
duality of descriptions and, hence, about the dual Meissner eifect-for loops instead of particles^^. 

Before describing the emerging picture in more detail, we note the following. Consider the expression 
for vorticity u) = k f dav{a)S{r — r{a)) from the point of view of reparametrization invariance. In 



particular, since we have a closed contour, we can always choose it as going from infinity to minus infinity 
(it is easy topologically to wrap it onto a closed contour of any size). For the function y = exp{a) we 
have evidently < y < oo when a varies from -oo to oo. This means that we can replace a by In y in the 
expression for the vorticity in order to obtain 



which in a nutshell is the same thing as a Dirac monopole, Ref.[70], with charge strength k, so that the 
vortices can be treated as Dirac monopoles. In Appendix C we provide some facts about Dirac monopoles 
in relation to vortices. According to Dirac [70] the string attached to such a monopole can either go to 
infinity (as in the present case) or to another monopole of equal and opposite strength. In our case this 
means only that when two hard spheres become hydrodynamically entangled, they cannot escape the 
linkage they formed. This is the (topological) essence of quark confinement in QCD known as monopole 
condensation^^. But we are not dealing with QCD in this work! How then we can talk about the QCD? 
The rationale for this was put forward first by Nambu, Ref.[71]^^. In his work he superimposed the G-L and 
Dirac monopole theories to demonstrate quark confinement for mesons (these arc made of just two quarks: 
quark and antiquark). For this qualitative picture to make sense, there should be some way of reducing 
QCD to G-L type theory. The feasibility of such an abelian reduction (projection) was investigated 
first by 't Hooft in Ref.[72]. Recent numerical studies have provided unmistakable evidence supporting 
the idea of quark confinement through monopole condensation, Ref.s [73,74]. Theoretical advancements 
made since the publication of 't Hooft's paper took place along two different (opposite) directions. In one 
direction, recently, Faddeev and Niemi found knot-like topological solitons using a Skyrme-typc! nonlinear 
sigma model and conjectured that such a model can correctly represent QCD in the low energy limit 
[75,76]. That this is indeed the case was established in a series of papers by Cho [77,78] and, more 
recently, by Kondo, Ref.[79]. In another direction, in view of the fact that, while macroscopically the 
Meissner effect is triggered by the effective mass of the vector field, microscopically this mass is generated 
by Cooper pairs [25], it makes sense to look at detection of the excited states of such Cooper-like pairs 
experimentally. The famous variational BCS treatment of superconductivity contains at its heart the gap 
equation responsible for the formation of Cooper pairs. The BCS treatment of superconductivity was 
substantially improved by Richardson, Ref.[80], who solved the microscopic model exactly. His model 
is known in literature as the Richardson model. Closely related to this model is a model proposed by 
Gaudin. It is also exactly solvable (by Bethe anstatz methods) [81]. The Gaudin model(s) describes 
various properties of one dimensional spin chains in the semiclassical limit. Energy spectra of the Gaudin 
and Richardson models are very similar. In particular, under some conditions they are equidistant, 
like those for bosonic string models. Recently, we were able to find new models associated with 
Veneziano amplitudes, e.g. see Ref.s [82,83], describing meson-meson scattering processes. In particular, 
we demonstrated that the Richardson- Gaudin spin chain model producing equidistant spectra can be 
obtained directly from Veneziano amplitudes. Since the Veneziano amplitudes describe extremely well 
the meson mass spectrum, and since we demonstrated that the Richardson-Gaudin model (originally used 
in superconductivity and nuclear physics) can be recovered from combinatorial and analytical properties 

^^It should be noted that in the case of usual superconductors one should distinguish between the constant magnetic fields 
penetrating superconductors and the fields made by vortices. In the case of colloidal suspensions it is also possible to create 
some steady velocity current and to consider velocity at a given point in the fluid as made of both steady and fluctuating 
parts. 

^^That is the Bose-Einstein-type condensation in view of results of Section 5. This is explained further in Appendix C 

^^We discuss his work briefly in Appendix C 

^*Also, for monopoles models discussed in Appendix C. 



C 








(6.31) 



— OO 
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of these amplitudes, this means that the Abehan reduction can be considered as confirmed (at least for 
mesons) not only numerically but also experimentally. 



6.6 Miscellaneous 



In Section 3.3. we demostrated that for colloidal suspensions it is sufficient to use only the Abelian version 

of the Chern-Simons theory for description of emerging entanglements. There could be other instances 
where such an Abelian treatment might fail. Examples of more sophisticated non- Abelian fluids were 
considered in several recent excellent reviews [84,85]. These papers might serve as points of departure for 
the treatment of more elaborate hydrodynamical problems involving non- Abelian entanglements. Finally, 
the force-free equation w = av which is used in our work, is known to possesss interesting new physical 
properties when, instead of treating a as a constant, one treats a as some function of the coordinates. 
Such treatment can be found in Ref. [86] and involves the use of conformal transformations and invariants 
recently considered in our work on the Yamabe problem, Ref. [87], and the Poincare' conjecture, Ref. [68]. 
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Appendix A. Some facts from the theory of Green's functions 
Consider an equation 



^ - ) $ = 0. (A.l) 



Such an equation can be written in the form of an integral equation as follows 



$(x,t)= J Go{x,t;x',t')<^>oix',t')dx'dt' (A.2) 

so that 

$(x,t ^ t') = ^'o(x,t')- (A.3) 
Under such conditions, the Green's function Go(x, t; x', t') must obey the following equation 

- Go(x, t; x', t') = 5{x - x')5{t - t') (A.4) 

provided that Go = for t < t' . In a, more complicated situation, when 

(^^^-H-V^ G(x, t; x', t') = 5(x - x')5{t - t') (A.5) 

we can write a formal solution for G in the form of the integral (Dyson's) equation 

G(x, t; x', t') = Go(x, t; x', t') + j Go(x, t; x', t')V{^, t')G{x', t'; x", t")dx'dt' (A.6) 

or, symbolically, G = Go + GqVG . In the case of Eq.(4.35) of the main text, we have to replace Eq.(A.5) 
by 

(J^-Hi-H2- Vi2^ G(xi,X2,i;xi,x^,i') = ,5(xi - x[)S{x2 - x'2)<5(i - t') (A.7) 
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and, accordingly, the Dyson type Eq.(A.6) is now replaced by the analogous equation in which now 
we must have Go(xi, X2, t; x'^, Xj, t') = Go(xi,t;x'j,i')G'o(x2,t;x2,t')- To check the correctness of such 
a decomposition we note that for x x' Eq.s (A.l) and (A. 4) coincide while for t ^ t' integration of 
Eq.(A.4) over a small domain around zero and taking into account that Go = for t < t' produces 
Go(x, t t';x',t') = (5(x — x'). Repeating those arguments for the two-particle Green's function and 
using Eq.(A.7) (with V12 = ) provides the needed proof of the decomposition of Go in the two-particle 
case. 

Define now formally the renormalized potential V via 

G = Go + GoVGo. (A.8) 
Then, by comparing this equation with the original Dyson's equation for G we obtain 

G - Go = GoVGo = GoVG = GoF(Go + GoVGo) (A.9) 
This allows us to write the integral equation for the effective potential V as 

V=V + VGoV. (A.IO) 



Appendix B. Dual treatment of the dynamics of colloidal supensions and hydrodynamic 
screening 

We begin by first considering screening. The path integral for the functional, Eq.(5.5), can be conve- 
niently rewritten as follows 



2we ^ 



(V-i— A)^ 



} 



= f/^M(VxAf+(^)'(V^-^An (B.l) 
upon substitution of the ansatz ip = —- ex.p{iip) into first line of Eq.(B.l). Such a substitution is 

ZTT 

consistent with the current defined in Eq.(4.50). 
Since V • A =0, we obtain 

( -'-^Af = mr + a) ^ 1^ A . V^. (B.2) 
Consider now the following path integral 

Z = l D{i;}eM-l[^y J dM(V^)'-§^A.V^)]. (B.3a) 

Since it is of a Gaussian-type, it can be straightforwardly calculated with the result 

Z = NeM-^\^A,). (B3b) 

Here N is some (normalization) constant. Using this result and Eq.(B.l) wc obtain the following final 
expression for the partition function for the vector A-field with account of constraints 

E = l D[A] exp{-^ I d^r{A^[-6,^V' - (1 - l)d^d.]A^ + e'A^{S^. - ^)A^}. (B.4a) 

This result is in complete accord with Eq.(4.51b) where for the mass m of the vector field A we obtained: 
m = e. The above derivation was made without the use of Higgs-type calculations, Ref [88] . Surely, it is 
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in accord with these calculations. We would like now to rewrite the obtained result in a somewhat formal 
(simplified) form as follows: 

E = l D[A]6{V-A)eM-^ j d'r[{V x Af + e^A^]}. (B.4b) 

This will be used below in such simplified form. To avoid extra notation, we also set , ^ = 1. This 

factor can be restored if needed. 

Now we are ready for the dual treatment, which can be done in several ways. For instance, following 
the logic of Dirac's paper [70 ], we replace S by 

E = J D[A]S{V ■ A) exp{-^ J dh[ (V x A) + v)^ + e^A^]} (B.5) 



a; 



where v = — = <f da-v{a)6(r — r(cr)). Next, we use the Hubbard-Stratonovich-type identity allowing us 



e 



c 



to make a linearization, e.g. 

ip{-^ j d^r[{ (V X A) + v)^} = J exp[-i J (i^r*^ + i J d^r{ (V x A) + v) • *] (B.6) 



exr 



Then, we take advantage of the fact that (V x A) • = (V x 'S') • A + V • (A x *).By ignoring surface 
terms this allows us to rewrite the above result as follows 

j exp[-^ J d^r^^ + i J d^r{ (V x A) + v) • *] 

= J exp[-^ j rf^r*^ + i j d^r{ (V x *) • A + v *)] (B.7) 

Using this result in Eq.(B.5) and using the Hubbard-Stratonovich transformation again we obtain: 

S = y D[A]SiV ■ A) exp{-i J d^rf (V x A) + v)^ + e^A^]} 

= J D[*](5(V • *) exp[-i j d^r[^^ + ^ x *)^] + i j d^^^ ■ v]. (B.8) 

Since exp(i J d^r^ • v) = exp(i f dav{a) ■ *(ct)) we can use this expresiion in Eq.(5.34) in order eventually 

c 

to arrive at the functional of G-L-type (analogous to Eq.(5.6) with obviously redefined constants). The 

vector field ^ is now massive. It is convenient to make a replacement: to make the functional 

for the "if field look exactly as in Eq.(B.5) (with v = 0). The above transformations provide a manifestly 
dual formulation of the colloidal suspension problem. These transformations can be made differently 
nevertheless. Such an alternative treatment is useful since the end result has relevance to string theory 
and to the problem of quark confinement in QCD as was first noticed by Nambu, Ref.[71]. This topic is 
discussed briefly in the next appendix. 



Appendix C Nambu string and colloidal suspensions: Some unusucd uses of Dirac 
monopoles. 
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We begin with Eq.(B.5) but this time we treat it differently. In particular, we have 



S = JD[A]6{V- 


■A) 


exp{ 


= Jd[A]S{V- 


■A) 


exp{ 


= Id[A]6{^- 


■A) 


exp{ 


= [d[A]S{V- 


■A) 


exp{ 



The exponent in Eq.(C.l) is useful for comparison with that given in Eq.(4.40). Such a comparison 
suggests that while the second (linking) term is essentially the same as in Eq.(4.40)^^, the first term in 
the exponent of Eq.(C.l) might be analogous to the "kinetic" string-like term in Eq.(4.40). This line of 
reasoning can be foimd in the paper by Nambu [71]. If one ignores quark masses as is usually done in 
string-theoretic literature, then Eq.(13) of Nambu's paper looks very much like our Eq.(C.l), provided 
that we identify the first term with the stringy Nambu-Goto term^^. To do so, we formally need to use 
the results of our Sections 5.5 and 6.2. This time, however, we have to allow for self-linking. Also, we have 
to take into account that for this case the energy and the helicity become the same (up to a constant). 
Thus, one can consider the helicity instead of energy. A very detailed treatment of helicity was made in 
the paper by Ricca and Moffatt, Rcf.[89], from which it follows that the helicity is ideally suited for the 
description of self-linking. In such a case we have to deal with closed curves of finite thickness. In fact, it 
is sufficient to have a closed tube instead of a closed infinitely thin curve. On such a tube one can perform 
the Dehn surgery by cutting a tube at some section, twisting the free ends through a relative angle 27rno, 
where no is some integer, and reconnecting the ends. This operation makes a self-linking proportional 
to riQ. If we agree that the Dehn twists are made only in increments of ±27r, we obtain the "spectrum" 
which is equidistant and, hence, string-like. 

This intuitive picture can be made more quantitative as follows. Taking into account Eq.(5.48), the 
kinetic term in the exponent of Eq.(C.l) can be tentatively written as follows 



' Ci Ci 

This expression suffers from two apparent deficiencies. First, while the second term in the exponent of 
Eq.(C.l) accounts for screening effects, Eq.(C.2) is written without such an accoimt. Second, since energy 
and helicity are proportional to each other and since the Dehn surgery can be made only for surfaces, 
Eq.(C.2) should be modified by replacing infinitely thin contours by tubes. To repair the first problem 
we follow the book by Pismen, Ref.[90], where on page 186 we find the following information. Consider 
our Eq.s(4.52a) or (4.52b) and take into account Eq.(4.15). Then, we can write 



V^A - e^A = -e dav{a)S{r - r{a)). (C.3) 



^2 A ^ _^ 

c 



The solution of the equation for vector potential A, Eq.(5.38), should now be modified to account for 
screening and boundary effects. The result for energy, Eq.(5.47), now will be changed accordingly so 
that the screening exponent will emerge in Eq.(C.l). To account for surface effects we recognize that the 
self-linking expression, Eq.(C.2) is reparametrization invariant. If, instead of infinitely thin contours we 
consider fiuctuating tubes, the reparametrization invariance should survive. The surface analog of the 
expression ^ dav{a)S{r — r(cr)) is given in Eq.(6.16b). By introducing the notation 
c 

da dr dr da 



■^^We have mentioned already that the screening is not affecting the topological nature of this term. 
26E.g. see Eq.(6.14). 
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the self-linking term can be brought into the following final form (for just one loop for brevity) 



which is just what Nambu obtained. He further demonstrated that such a term can be transformed 
into — m / dadr^/^ (e.g. see our Eq.(6.14)) with the constant m (the string tension) being related to 
coupling constant(s) of the theory. Since in the limit of infinitely thin tubes results just obtained match 
those discussed in our Section 5.5, wc would like to take advantage of this observation. In Section 5.5 we 
considered fully flexible (Brownian) loops. From the theory of polymer solutions it is known that such 
loops can be made of the so called semiflexible polymers whose rigidity is rather weak. Following our 
work, Ref.[91], the path integrals describing semiflexible polymer chains are given by 

I = j D[n{T)]eM-SHT)]) (C.6) 

with action S'[u(r)] given by 

TV 2 ^ 

SHt)\ = f /'^^ (^) + / dTX{T){u\T) - 1). (C.7) 



The rigidity constant is k. For brevity it will be put equal to one. The Lagrange multiplier A takes care 
of the fact that the "motion" is taking place on the surface of a 2-sphere. Minimization of the action S 
produces 

= Au (C.8) 

with Lagrange multiplyer being determined by the constraint = thus producing instead of Eq.(C.8) 
the following result: 

u = -(u-u)u, (C.9) 

where u =^u, etc. In view of the results of this subsection, consider now an immediate extension of the 
obtained results known as the Neumann modeP^ 

u+Gu=Au, u2 = l (C.IO) 

for some matrix G which always can be brought to the diagonal form. By analogy with Eq.(6.8), we can 
rewrite Eq.(C.lO) in the following equivalent form 

ux[u + Gu]=0 (C.ll) 

since u x Au =0 . The above equation is just a special case of the Landau-Lifshitz (L-L) equation describing 
dynamics of Heisenberg (anti)ferromagnets. In one space and one time dimension the L-L equation reads 

|u = {ux[u + Gu]}, (C.12) 

where now u =^u, etc. In Sections 6.3 and 6.5 we mentioned already that L-L equation describes the 
dynamics of vortex filaments in fluids, plasmas, etc and is also obtainable from the Lund-Regge theory. 
Following Veselov.[93], consider a special solution of the L-L equation obtained by inserting the ansatz 
\i{x,t) = \i{x — iOt) into Eq.(C.12). Such a substitution produces: 

-z6'u={ux [u-HGu]},u^ = L (C.13) 



■^'^Some useful details related to Neumann's model can be found in our work, Ref.[92]. 
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Taking a vector product of both sides of this equation produces 

u + Gu =\u+ie[u X u] . (C.14) 

This equation describes the classical motion of a charged particle in the presence of a Dirac monopole. 
At the quantum level such a problem was studied in detail by Dunne, Rcf. [94], who demonstrated that 
in the limit 9 oo the monopole spectrum is equidistant. This result is compatible with the result of 
Ricca and Moffat [89] , and explains the role of monopoles in quark confinement (in view of results of our 
Section 5.5). Furthermore it corroborates the results of our recent work, Rcf. [83], briefly mentioned in 
Section 6.5., where the spectrum of thef-d Heisenberg XXX spin chain was recovered directly from the 
combinatorics of scattering data supplied by uses of Veneziano amplitudes in scattering experiments. 
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